ANISOTROPIC CONTACT PROCESS ON 
HOMOGENEOUS TREES 

Irene Hueter 
Abstract 

The existence of a weak survival region is established for the anisotropic symmetric 
contact process on a homogeneous tree of degree 2d > 4 : For parameter values 
in a certain connected region of positive Lebesgue measure, the population survives 
forever with positive probability but ultimately vacates every finite subset of the tree 
with probability one. In this phase, infection trails must converge to the geometric 
boundary of the tree. The random subset A of the boundary consisting of all ends of 
the tree in which the infection survives, called the limit set of the process, is shown to 
have Hausdorff dimension no larger than one half the Hausdorff dimension of the entire 
geometric boundary. In addition, there is strict inequality at the transition between 
weak and strong survival except when the contact process is isotropic. It is further 
shown that in all cases there is a distinguished probability measure /i, supported by fi, 
such that the Hausdorff dimension of A n , where fi^ is the set of ^-generic points of 
fl, converges to one half the Hausdorff dimension of at the phase separation points. 
Exact formulae for the Hausdorff dimensions of A and A n fi^ are obtained. We also 
prove that the contact process at the transition between extinction and weak survival 
does not survive. The method developed shows that the contact process at the phase 
transition to strong survival survives weakly for d > 2. 

1 Introduction 

1.1 Background: Transition from Weak to Strong Survival 

This paper considers the anisotropic contact process on an infinite homogeneous tree when 
in weak survival or at the transition to strong survival. The process was introduced in 
1 23 1 and pursued intensely in the isotropic briefly surveyed in ||2^ (see |l^ 



for the contact process in general). In fact, the contact process is a stochastic growth 
process which, along with branching random walks and percolation processes in spaces with 
hyperbolic geometries, exhibits an intermediate phase not present in the corresponding 
processes living in spaces with Euclidean geometry. This is the weak survival region, in 
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which the "population" survives forever with positive probabiUty but, with probabiUty one, 
eventually vacates every compact subset of the ambient space. Weak survival is known to 
arise for the isotropic contact process on a homogeneous tree |23, 27], for the anisotropic 
symmetric branching random walk on a homogeneous tree |^], branching Brownian motion 
in the Poincare plane |14], and site percolation on a co-compact Fuchsian group [12|. The 
transition from weak to strong survival still eludes complete understanding for many particle 
systems. 

This paper establishes the existence of a weak survival region and analyzes its features 
for the anisotropic symmetric contact process on an infinite homogeneous tree of even de- 
gree > 4. For ease of exposition, we shall restrict our attention to the contact process on 
the infinite homogeneous tree of degree 2d for d > 1. Applying our technique to the 
anisotropic nearest neighbour contact process with some symmetry assumption on homo- 
geneous trees of odd degree is more subtle (see Remark (4) at the end of the Introduction) 
and leads to slightly different algebra. 

The anisotropic case poses some genuine difficulties not present in the isotropic case. The 
approach pursued here is completely different from the existence proofs for weak survival as 
previously given for the isotropic process. For instance, the proof in this paper does not use 
the weight function and the fact that the contact process becomes extinct at the transition 
to survival (even though the latter result will be proved subsequently). Our main interest 
is to study the limit set A, defined to be the subset of the geometric boundary (the set 
of ends) of the tree in which the infection survives, when the contact process is in weak 
survival and to explore the behaviour of the transition between weak and strong survival. 
As our explicit formulae show, the Hausdorff dimension 6h (A) of the limit set never exceeds 
i the Hausdorff dimension 6h{^) of the geometric boundary, and importantly, the equality 
(A) = ^6h{^) is valid exactly at the transition between both survival regions if and only 
if the contact process is isotropic. This confirms a conjecture, raised in |15|, that there 
is equality in the isotropic case. Moreover, we will prove that, at the transition between 
extinction and weak survival, the contact process on 7^^ does not survive for d > 1. Our 
method shows that at the phase transition to strong survival, the contact process on 
survives weakly for d > 1. 

Additionally, we will investigate (a) the exponential rate r] of decay in time t of the prob- 
ability that the initial infected site is infected at time t, (b) the exponential rate of growth 
in space-time on the event of survival, (c) the distances of the nearest and furthest infected 
vertices from the root vertex at time t on the event of survival, and (d) the distribution 
of the limiting points in weak survival. Key ingredients to our analysis are shift-invariant 
probability measures supported by O and their associated transition matrices related to a 
stationary one-step Markov chain on some finite set of generators. These probability mea- 
sures arise from "normalizing" the infection probabilities at each large distance from the 
root vertex as the infection is moving off to the boundary of the tree. 

This paper is an expanded version of talks that I presented at an AMS Special Session 
in Gainesville, Florida (March 12-13, 1999), at a Probability Meeting at Colorado Springs 
(May 28-30, 1999), at the World Congress of the Bernoulh Society and IMS Meeting in 
Guanajuato, Mexico (May 15-20, 2000), at an AMS Special Session in New Orleans (January 
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13, 2001), at a DIMACS/DIMATIA Workshop at Rutgers University, New Jersey (March 
19-21, 2001), and in Colloquia or Seminars at the New College of the University of South 
Florida, Sarasota (October 1, 1998), at the University of Wisconsin at Milwaukee (November 
13, 1998), at Duke University, Durham, North Carolina (September 15, 2000), and at the 
University of Berne, Switzerland (April 23, 2001). 

1.2 Anisotropic Contact Process 

We shall restrict our attention to an anisotropic, symmetric contact process on a homoge- 
neous tree of even degree 2d {d> 1). An anisotropic contact process on the tree T = 
is a continuous time Markov process At on the set of finite subsets of (the vertex set of) T 
that evolves as follows. Infected sites (members of At) recover at rate 1 and upon recovery 
are removed from At- Healthy sites (members of Af) become infected at a rate that equals 
the sum of the infection rates attached to the edges leading to infected nearest neighbours 
and upon infection are added to At- Under the default probability measure P, the initial 
state ^0 is the singleton set {1} (where 1 is a distinguished element of T, called the root). 
Each vertex a; of T has exactly 2d neighbours. The tree T is homogeneous in that for any 
two vertices x and y there is an isometry that maps x to y. Associated with each of the 
2c? emanating edges is an infection rate. The symmetry assumption guarantees that each 
infection rate be used twice for the set of emanating edges of each vertex. In particular, the 
same infection rate is attached to an edge when the "infection crosses the edge" forwards 
and backwards. Hence, in the notion explained in Section 1.3 below, the set of infection 
rates is Am, . . . , Aa^, A^-i, . . . , A^-i, and for each of the d letters j, we assume \j = A^-i. 
If the infection rates are all equal, then the contact process is called isotropic. Note that, 
alternatively, in considering the question of existence of weak survival for the non-isotropic 
contact process one may vary the recovery rates instead of the infection rates to define an 
anisotropic process. 

1.3 The Tree as a Cayley Graph 

Representing the tree T = as the Cayley graph Q = on d generators proves its worth 
to understand the anisotropic contact process on the tree T. Let .4+ = {ai, 02, . . . , a^} be 
a set of d letters, let A- = {a]"^,a^''^, . . . ,a^^} be the set of formal inverses of the letters 
in A+, and set A = A+ U A-- The free group Q with generators A+ is the set of finite 
reduced words from the alphabet A (a word is reduced if no letter a G A is adjacent to its 
inverse) , where multiplication is concatenation followed by reduction and the group identity 
1 is the empty word. There is a natural bijection between G and the set of vertices of T, 
in which g,h e Q are mapped to adjacent vertices of T if and only if gh~^ G A. In other 
words, vertices are uniquely represented by finite reduced words from A. In the subsequent 
discussion, we shall not be careful to distinguish between vertices of T and the words (or 
group elements) representing them, and we shall refer to G as the vertex set of T. For any 
vertex z, denote by \z\ the length of its representative word. Note that \z\ as well is the 
distance from vertex z to vertex 1 in the graph T. For every integer n > 0, let Qn denote 
the set of all vertices x e Q aX distance n from the root vertex (i.e. |x| = n). 
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In a canonical way, the bijection between T and Q induces a bijection between the natural 
boundary of T and the geometric boundary fl, the set of semi-infinite reduced words from 
the alphabet A. A geodesic in T is a finite or semi-infinite sequence of distinct vertices 
vi,V2, ■ ■ ■ such that for every i > 1, the vertices Vi and Vj+i are nearest neighbours. An end 
of T is an equivalence class of semi-infinite geodesies, two geodesies being equivalent if and 
only if the sets of vertices through which they pass differ in at most finitely many vertices. 
If a; = X1X2 ■ ■ ■ £ then oj corresponds to the end of T represented by the semi-infinite 
geodesic that passes through the vertices 1, xi, xi2;2, . . . in succession. For each real number 
a £ (0, 1), there is a natural metric da on 0,, defined by 

da{uj,u;')=a''^^'^'\ (1.1) 

where N{uj,uj') is the largest integer n such that the sequences uj and uj' agree in entries 
1,2, ... ,n. For any choice of a, the corresponding topology on 17 is the topology of coordi- 
natewise convergence. For any vertex z of T, define T{z) to be the set of vertices v such that 
the geodesic segment from 1 to u passes through z, equivalently, such that the unique word 
representing z is a prefix of the word representing v. Similarly, define f2(z) to be the set of 
infinite reduced words iv = xiX2 ■ ■ ■ such that, for some finite n, the word z is represented 
by the word X1X2 • • • Xn- Observe that, for every integer n > 1, the set {^l{z) : \z\ = n} is a 
finite open cover of the geometric boundary 0,. Finally, define S to be the set of all doubly 
infinite reduced words ^ = {xn)^=-oQ from A. 



1.4 Anisotropic Contact Process on Q 

The symmetric contact process on a homogeneous tree of degree 3 or larger distinguishes 
itself from the symmetric process on the integer lattice Z*^ in that there are two different 



survival regions 19, 27]. More precisely, on T2d for d > 1 (but not for d = 1), there is a 
partition of the parameter space [0, oo)'* = R'J_ 9 A = {{Xa}aeA+) into three regions 7^i,7^2 
and TZs such that 

(a) if A € TZi, then = eventually, with probability 1, 

(b) if A G 7^2, then P{\At\ ^ 00} > 0, but Vx G T, P{x £ At for arbitrarily large t} = 0, 

(c) if A G TZs, then with positive probability \At\ 00 and, for all x G T, for arbitrarily 
large values of t, x £ At. 

On TZi, the contact process is called subcritical, on TZ2, weakly supercritical, and on TZ^, 
strongly supercritical. The main results of this paper concern the weak survival regime 

AG7^2, (1.2) 

whose existence needs to be established first. Additionally, our findings will shed some light 
onto the nature of the boundary of 7^2, that is, the phase transitions. 

We begin with introducing some terminology for the anisotropic contact process with 
infection rates {Xa}a£A+ ^-iid recovery rate 1, in order to state the main results. Thus, for 
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every x £ Q, vertex x attempts to infect vertex xi, i £ A, at rate Aj (note that vertex xi 
may already be infected). Recall that 1 denotes the root vertex. Define 

= r?(A) = lim (P{1 e At})^/^ = sup(P{l £ ^J)^/* < 1. (1.3) 

The limit exists because by the Markov property (see Section 2.1) and monotonicity prop- 
erties of the contact process, P{1 G At}P{l £ As} < P{1 £ At+s} for all real s,t > 0. 
Thus, an easy subadditivity argument applies. It follows as well that P{1 £ At} < ?/(A)* 
for alH > 0. It is obvious that the function r]{-) is nondecreasing in each Xj. For any vertex 
X £ T, define 

Ux = UxiX) = P{x £ At for some t > 0}. (1.4) 

The strong Markov property together with the monotonicity properties of the contact pro- 
cess and the homogeneity of the process at each vertex implies that u^y > u^Uy for each 
x,y £ T such that \xy\ = \x\ + \y\ (no reduction occurs when x and y are concatenated). 
If the contact process is weakly supercritical, then Ux < 1 for every x 7^ 1 in T because 
Ux = 1 would imply that, with probability one, the root be reinfected at indefinitely large 
times, thus, the process would be strongly supercritical. A subadditivity argument shows 
that, for every a £ A and each vertex x = aa . . . a £ Gn, the limit 

, hm UxiX)'^'' = Pa = Pa{X) (1.5) 

exists and that Ux{X) < /?a(A)" for all n > 0. Moreover, for every integer k and x £ 
Qk, by subadditivity (more precisely, by supermultiplicativity) , for each periodic sequence 
yn = XX . . . x £ Qnk ■, the limit 

Jim n,„(A)V-=/3,(A) (1.6) 

exists for every A and Uy^{X) < /^^..(A)" for every n > 0. Clearly, the functions Pxi') are 
nondecreasing in each infection parameter Aj. We will discuss the strict monotonicity prop- 
erties of the PxW, vWi and other functions (Sections 2.5 and 4.1) and their continuity 
properties (Section 5). 



1.5 Main Results: Weak Survival Region and Limit Set 



There has been a wealth of results on isotropic contact processes on trees, among them 
|2|, I, m, 13, |18|, |0|, H, [l|, m. In contrast to the contact ] process on Euclidean lattices 
the isotropic contact process on homogeneous trees was shown to exhibit two essentially 
different survival phases, the weak one being a novelty |23, |l^]. One might suspect the 
anisotropic contact process to as well have an intermediate region between extinction and 
strong survival. A result of this paper confirms this previously conjectured behaviour for 
the symmetric anisotropic contact process on T2d- 

Further detailing the boundary of the weak survival phase involves functions < bi{X) < 1 
{i £ A), nondecreasing in each argument Aj, which unfortunately do not come with a short 
description but are defined in a number of steps and merely are intermediate tools in our 
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study (as opposed to being of interest on their own). However, to characterize the hmit 
set at the weak/strong survival transition, the 6j(A) provide a crucial link to an analytic 
tool available for the symmetric anisotropic random walk The 6j(A)^ arise as entries of 
a certain 2d x 2d Perron-Frobenius matrix with a lead eigenvalue that coincides with the 
exponential of the expectation of a potential function relative to a Gibbs measure. The full 
description is deferred to Sections 3 and 4. Importantly for the next result, each < 6i < 1. 

Theorem 1 The weak survival region IZ2 is nonempty unless d = 1 and enjoys the following 
properties: 



(a) There are functions bi{X), defined in such that the boundary TZi H 7^2 consists 

of all A for which the 6j(A) satisfy 



and such that the boundary 7^2 H TZ^ consists of all A so that 

sr- h{X? _ 

(b) Every line in the interior of the first quadrant in R'^ that passes through the origin 
has an intersection with IZ2 that is a line segment. 

(c) The region IZ2 has positive d-dimensional Lebesgue measure. 

(d) The region IZ2 is connected and is a symmetric region in the d parameters Aj. 

The critical contact process behaves as follows. 
Theorem 2 For A € IZ2 H 7^3, the contact process survives weakly. 
For a proof of this result, see Corollary ^ 

Theorem 3 T^-i H 1^2 C T^i, that is, for A G n IZ2, the contact process on 7^^ for d> 1 

almost surely becomes extinct. 



In the isotropic case, the result is in for degree 2d > A and in for degree 3. 
Moreover, the critical contact process on Z and in all other Euclidean lattices dies out Q. 
Theorems |l] through |3| will be proven in Section 5.6. For any set B C R*^, let B'^ denote its 
complement in R"^. 

Theorem 4 For all d > 1, we have r]{X) < 1 for X G int(lZi U IZ2), and, r]{X) = 1 for 
XGV^nns. 

In the weak survival region TZ2, the "population" eventually vacates every finite subset 
of vertices of the tree with probability one. Therefore, the population has a well-defined 
limit set. Define the limit set A of the contact process on T to be the (random) set of 
uj = X1X2 . . . S such that each vertex X1X2 ■ ■ ■ x^ of w is infected at some time. It is 
easily seen that, if the contact process is supercritical, then on the event of survival, A is 
nonempty and compact (relative to any of the metrics da)- 
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Theorem 5 For A G 7^2 1 almost surely on the event of survival, the Hausdorff dimension 
6{\) of A (relative to the metric da, defined in ( [i. is given by 

m--'^, (1.7) 

log a 

where 6{X) is the leading eigenvalue of some Perron-Frobenius matrix, described in Section 4, 
and is the unique positive number such that 

where the 6j(A) are defined in For X £ Tl2, the functions 6{X) and6{X) are continuous 

functions in each of the variables Xj and are strictly increasing along "directions of increase" 
for X (see Section 2.5). Furthermore, if 5h{^) denotes the Hausdorff dimension ofQ, 

S{X)<^5Hm, (1.9) 

with equality holding if and only if the underlying contact process is isotropic and X G 7^2n7^3 
(at the transition between weak and strong survival). 



The equality in ( |1.9| ) confirms a conjecture, first raised in concerning the behaviour 
of the isotropic contact process at the transition between weak and strong survival. Any 
function in A that is strictly increasing along "directions of increase" for A necessarily is 
strictly increasing if every component of A is increased. The following properties will be 
valuable. 

Theorem 6 For each i,j G A, the functions 

Xj ^ 6i(A), A, ^A(A), 

A,- ^ rj{X), A, ^ 0(A) 

are continuous for each X G int(TZi U 7^2)5 are left- continuous for X TZ^ such that each 
Xk > 0. Moreover, on int(TZi U 7^2)1 each of these functions is strictly increasing along 
"directions of increase " for X as defined in Section 2. 5. 

Now we turn back to the limit set of the contact process. The more intriguing part of 
our studies revolves around the behaviour of the limit set A when A takes a critical value in 
7^2 n 7^.3. In order to understand the fine structure of A, we partition the geometric boundary 
Q into measure classes fi^ and ask about the size of the intersection of each equivalence 
class with the limit set. Recall that Q is the set of semi-infinite reduced words from A. Let 
o" : — > be the one-sided forward shift operator on 0,, that is, 

a{xiX2 ■ ■ ■) = X2X3 

For any ergodic, cr-invariant probability measure fi on J7, define to be the subset of 0, 
consisting of all a; G such that for every continuous real- valued function / : $7 — > R, 

. n 

hm - ^ fia'u;) = / fdfi. (1.10) 
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Birkhoff's ergodic theorem implies that /u(r2^) = 1 since the space of continuous functions 
on 17 is separable in the sup norm topology. Moreover, if /x and v are distinct ergodic 
probability measures, then Q,^r\VLy = 0. 

For any ergodic, a-invariant probability measure /i on 0, let h{ii) denote the Kolmogorov- 
Sinai entropy of the measure-preserving system (r2,/x,cr) (for the definition, see e.g. |28|, 
Chapter 4). Define the function (px : Q. ^ Yihy 

ipx{xiX2...) = \ogK^{\), (1.11) 

where the functions 6j(A) are defined in ([4.1]). 

Theorem 7 Let A G 7^2 o-nd let fi be any ergodic, a-invariant probability measure on Q. If 
h{ii) < — J fxdii, then almost surely, A n $7^ = 0. // /i(/i) > — J fxdfi, then almost surely 
on the event of survival, the set An 17^ has Hausdorff dimension 6{X; fi) (relative to the 
metric da) 

s^x.^) = -M±M^, (1.12) 

log a 

This Hausdorff dimension satisfies the inequality 



6{X;f,)<-5Hin^), (1.13) 



where equality holds in for one and only one ergodic probability measure fj,^ and only 

when A G 7^2 n 7^3. 

The second display is the more curious part of Theorem 0, the existence of a shift- 
invariant probability measure /U* on O for which equality holds in (|1.13| ), and provides a 
more vigorous instance of equality (|1.9D , its analogue in the isotropic case. In fact, the 
transition from weak to strong survival happens precisely when, for some O^, the set A fills 
a subset of half the Hausdorff dimension of O^. The distinguished probability measure /i* 
is defined as follows. 

Define the backscatter matrix M2 = M2(A) to be the 2dx 2d matrix, indexed by elements 
of A, whose entries are given by 

(M2(A)),, = 6,(A)2 iij^i-\ (1.14) 
= ifj = i-i. 

We will see that, for A G TZ2 such that each A^ > 0, this is an irreducible nonnegative 
matrix, thus, the Perron-Frobenius theorem applies. Proposition 19 below will show that 
the lead (Perron-Frobenius) eigenvalue of M2 is 1 for every critical value A G TZ2 n TZ^. It 
follows that, if v is the (positive) right eigenvector, then for every A G 7^2 H TZ^, 

P2{hj) = — (1.15 

are the entries of an irreducible stochastic matrix P2. The probability measure is the 
unique probability measure on O such that the induced coordinate process is the stationary 
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Markov chain with transition probabiHty matrix P2. In fact, it is only possible that such an 
invariant measure exists for each A at the phase transition 7^2 H TZ^ due to the fact that the 
critical contact process survives weakly. Observe that in a similar fashion, for each A € 7^2 
and each positive real p above some (critical) value, if Mp represents the matrix that results 
when each power 2 is replaced by power p in ( |1.14| ), then 

P,(i,i) = |^^(i-*.(r')) 

8{p; X)Vi 

denotes the transition probabilities of a stationary one-step Markov chain, where v is the 
right eigenvector associated with the lead eigenvalue 9{p; •) of the Perron-Frobenius matrix 
Mp and 5 (') denotes the Kronecker delta function. For each A G 7^2 and each such p, there 
exists a unique shift-invariant probability measure supported by Q (see Section 3 for more 
details) . 



Moreover, it will be demonstrated (Theorem 10) that, if 0^10^2 • • • is a limit point in 



and if pn denotes the distribution under P of the process tOn, ^n+i, ■ ■ ■ , then for every n > 1, 
the measure pn is absolutely continuous with respect to a Gibbs state p^ and pn Pip as 
n — > 00. The limiting distribution will be identified. The stationary process induced by p,p 
is isomorphic to a Bernoulli shift These measures decay exponentially in the distance 
from the root vertex and are spherically symmetric if and only if the contact process is 
isotropic (see also Proposition ^ in Section 3.5 as well as the remarks in the paragraph 
thereafter) . 

Our final result offers insight into the dispersal behaviour of the infection in space-time. 
The symmetric anisotropic contact process moves at linear distance with time as does the 
isotropic contact process [^]. Following the same notation, we define rt and Rt to be the 
smallest and largest distances among the infected sites x £ At and Nn{ns) to be the number 
of vertices x £ Ans at distance n from the root that are infected at time ns. In the time- 
dependent case, let exp{$i_<j.;^} be the analogous function to ^(A), defined in ( |1.8D , at time 
scale s, thus, exp{$i g.;^} is the unique positive number such that 

y — _M^1^ — = 1, (1.16) 

exp{$i,,;A} + 6i,.(A) 

where the hi^s are the time-dependent equivalents of the functions 6j, encountered in Theorem || 
(for a fuller description, see ( |3.37| ) and ( |3.38| )). In Q, the function V{s) = — log(2(i— 1) 



was called the growth profile of the contact process. 

Theorem 8 Let d > 1 and A in 1Z2- Then there exist smallest and largest solutions < 
si < S2 < 00 of '&i,s;A = 0. Almost surely on the event of survival, 

lim rt/t = l/s2, (1.17) 

I— >oo 

Rt/t = l/si. (1.18) 
Moreover, for each s > such that <l>i,s;A > 0, almost surely on the event of survival, 

lim - log Nn{ns) = ¥i,,.a. (1.19) 

n^oo 77, 
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See Section 10 for the proof. 



Remarks. 

(1) The phenomenon that the lead eigenvalue of M2 converges to 1 at the phase transition 
will be shown to be responsible for both convergences of the Hausdorff dimension in 
the lower bounds to ^ of some Hausdorff dimension in ( |1.9D and ( |1.13| ), and, offers a 
useful characterization of the phase transition. 

(2) Analogous results can be proved for a contact process with infection rates between 
sites within finite distance (but not exclusively nearest neighbours). 

(3) It is worthwhile noticing that the formulae for the Hausdorff dimensions resemble the 
ones for the branching random walks on T2d-, with the average asymptotic infection 
probabilities hi corresponding to the generating functions Fj. Surprisingly, for the 
isotropic models, the dimensions of the limit sets coincide at the phase separations. 
For each anisotropic contact process at the phase transition, there is an anisotropic 
branching random walk with the same Hausdorff dimension of the limit set. 

(4) How to define a suitable symmetry assumption for the contact process on a homoge- 
neous trees of odd degree is a bit less obvious. Suppose that, at any vertex, associated 
with each of the d emanating edges that point away from the root vertex be an in- 
fection rate Xi, i G A+, and associated with the emanating edge that points back 
towards the root be an infection rate A^-i for k £ A+. Thus, two infection rates 
and Afc-i are attached to an edge, one being in use when the infection "crosses the 
edge forwards" and the other being in use when the infection "crosses the edge back- 
wards" . Suppose that the contact process on the tree T^+i of degree d+1 > 2 satisfy 
the symmetry assumption Afc = Xk-^ for each k G 

Question 1: Which of our results continue to hold for this model ? 

We remark that this setup gives rise to nonhomogeneous vertices, more precisely, the 
set of infection rates is not the same at different vertices. The method developed in 
this paper relies on the geometric decay (see Lemma ^), which in turn invokes the 
homogeneity of the vertices. 

Preliminary results that we pursue in subsequent papers indicate that if A^^-i = A* 
for all k G A+ (isotropic backtracking), then weak survival exists and, for A G 7^2, 
all of Theorems ^ and |^ hold, where 9{X) = J2ieA+ '^iW^ with the possible exceptions 
of the second parts in (|1.9| ) and ( 1.13| ). Moreover, the discontinuity set consists of all 



A which satisfy the equation J2i£A+ ^iW'^ Aj/A* = 1. 

Question 2: Does the contact process on T^^i for d + 1 > 2 with A* < maxjg_4^ Aj 
always exhibit weak survival, more specifically, is the transition from TZi to TZ^ pre- 
cluded ? 

The rest of the paper is organized as follows. Section 2 discusses the percolation struc- 
ture, the Markov property, and the strict monotonicity of the (3x, introduces the embedded 
labelled Galton- Watson trees, and outlines the strategy of proof (Section 2.5) for the ex- 
istence of weak survival. In Section 3, we give some background on the Thermodynamic 
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Formalism, describe the potential functions and Gibbs states, and prove when ij < 1. Section 
4 defines the first-passage and backscatter matrices and analyzes the strict monotonicity 
of their lead eigenvalues. Section 5 is concerned with the continuity features of the lead 
eigenvalues, derives a characterizing equation at criticality, establishes weak survival, and 
includes the proofs of Theorems ^, ^ and Theorems [l| through |3| on the critical contact pro- 
cess. Section 5 also includes a result on the distribution of the limit points of the contact 
process (Theorem |l0| ) . The upper and lower bounds for the Hausdorff dimensions are taken 
care of in Sections 6 and 7. Section 8 verifies inequalities (^) and ( 1.13|) as an appeal 
to the Gibbs Variational Principle, which will complete the proofs of Theorems ^ and |^. 
Section 9 reviews the isotropic case and Section 10 addresses the dispersal of the contact 
process over space-time. 



2 Basic Properties and Embedded Galton- Watson Trees 

First we describe the strong Markov property and the visualization of the contact process 
as a percolation structure, upon which the verification of the strict monotonicity of the 
functions Ux along certain directions in the space of infection parameters is based. This 
feature together with the continuity properties will be essential to our proof of Theorem |l|. 
Furthermore, we will describe the Galton- Watson trees embedded in the set of vertices ever 
to be infected. 



2.1 Percolation Structure and Strong Markov Property 

Percolation structure. The contact process may be constructed via the usual percolation 
structure on T x (0, oo), that IS , ctS Si system of independent Poisson processes attached to 
vertices and ordered pairs of neighbouring vertices. For each vertex x £ T, the Poisson 
process attached to x has rate 1, and determines the recovery times, specifically, at every 
occurrence time, site x recovers if it is infected. For each ordered pair (x,xi), i £ A, of 
neighbouring vertices, the Poisson process attached to {x,xi) has rate Aj, the occurrence 
times being precisely those times when an infection at x may jump to xi. Occurrences in 
these Poisson processes are marked on a system of directed rays {x} x [0, oo) connected to 
the vertices x of T, so that (A) at each occurrence time t of the Poisson process attached 
to {x,xi) an infection arrow is drawn from {x,t) to {xi,t), and (B) at each occurrence time 
t of the Poisson process attached to x a recovery mark * is attached to {x,t). There are 
no simultaneous occurrences of infection arrows and/or recovery marks in the percolation 
structure. At time t, the contact process now consists of all those vertices y for which there 
is a (directed) path through the percolation structure, the system of rays and arrows just 
described, that begins at the root vertex 1, ends at {y,t), and does not pass through any 
recovery marks *. An infection trail is a connected path in the percolation structure that 
does not pass through any recovery marks. 

Strong Markov Property. An important property that the contact process enjoys is 



the strong Markov property (see also ||l5|]. Section 2.2). Let G,Fi,F2, ■ ■ ■ ,Fk be pairwise 



nonoverlapping parts of the tree T, and let Si, S2, ■ ■ ■ , Sk be stopping times determined 
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by the percolation structure over G. Then conditional on the percolation structure over G, 
the post-Sj portions of the percolation structures over the sets Fi are independent, and for 
each i, the post-5j percolation structure over Fi has the same distribution as the entire 
percolation structure over Fj. 



2.2 Downward Infection Trails 

Fix a generator a £ A, and consider the subtree T* = T— T(a~^) of T, each of which vertices 
except the root vertex 1 is represented by a word x G G beginning with a letter xi / a~^. 
We arrange the tree T* in levels Co,Ci,C2, ■ ■ ■ , where x € if and only if |x| = n. There 
are {2d — 1)" vertices at the nth level Cn- Moreover, for every n > 1, define C*^ to be the 
subset of Cn containing those vertices whose word representation x = xiX2 • • • x„ terminates 
in the letter Xn = a. Observe that for any vertex x G £* , the set of nearest neighbours of x 
in Cn+i is {xy : y G Ci}, thus, for x £ £*, the tree T{x) is the left translate of the tree T* 
by the group element x. 

Let X be a vertex in level Cn and let y be a vertex contained in T{x). In other words, 
the word x is a prefix of the word y and y must lie in a level Cn+m at larger distance than x. 
Define a downward infection trail from x to y to be an infection trail that begins at x, does 
not exit T{x), and first reaches Cn+m at y, where it terminates. For every vertex x £ Q, 
define 

"Dx = {3 downward infection trail {root} x beginning at t = 0}. 
Thus, if we write 

Wx = P{Vx}, (2.1) 

then as an appeal to the monotonicity and the strong Markov property of the contact 
process, we have Wxy > WxWy for every x,y G Q so that \xy\ = \x\ + \y\. It is obvious 
that Wx < Ux- For each periodic sequence yn = xx . . . x S Gnk with x £ Gk, the limit 
limn^oo{wy„{X)y^'^ exists for every A because of the homogeneity of the process at each 
vertex, and in fact, by the same arguments as presented for the isotropic contact process in 



|15|, 



lim(«;,„(A))V« = /3,.(A), (2.2) 



where /3x(A) was defined in (|l.6D . Moreover, Wy^{X) < l3x{X)"' for every n > 0. The non- 



periodic analogue of (2.2) in the anisotropic case will be derived in Section 5.2. 



2.3 Labelled Galton- Watson Processes and Trees 

Every Galton- Watson chain has its associated genealogical tree r. This Galton- Watson tree 
may be described as follows (see also [@]). Vertices of r are arranged in levels Vq, Vi, V2, . . . . 
The vertices of level Vn represent the individuals of the nth generation of the corresponding 
Galton- Watson chain. Edges of the tree connect children and their parent, thus, there are 
edges only between vertices of successive levels. The limit set Aqw of the tree is the set of 
ends, i.e. the set of all infinite paths that start at the root Vq and visit each of the levels 
Vn exactly once. If the Galton- Watson chain dies out, then the limit set is empty. For each 
a G (0, 1), the da-distance between two ends 7, 7' is defined to be a", where n = n(7, 7') is 
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the last level Vn where the paths touch. The set Aqw is a compact metric space with this 
metric da- 

For any finite set B, a labelled Galton-Watson process with label space B is determined 
by a probability distribution Q on the set of 2^ subsets of B. Each individual C, regardless of 
its type, produces a random set of offspring, with distribution Q, and the offspring sets 
of different individuals are conditionally independent, as in an unlabelled Galton-Watson 
process. Thus, a labelled Galton-Watson process is a multitype Galton-Watson process in 
which (i) all types have the same offspring distribution, and (ii) the offspring distribution 
is constrained to allow at most one individual of each label. If Zn denotes the cardinality 
of the nth generation, then the sequence {Z„}n>o is an ordinary Galton-Watson process. 
Moreover, if the expected cardinality of a random set chosen according to the distribution 
Q is bigger than 1, then {Z„}„>o is supercritical. 

A labelled Galton-Watson tree r is the genealogical tree associated with a labelled 
Galton-Watson process so that the labels of the corresponding individuals are assigned 
to the vertices of r. If the underlying Galton-Watson process is supercritical, then with 
positive probability the tree r is infinite. On this event, each end of r will be naturally 
identified with a unique semi-infinite sequence to = xiX2 . . . , where the end crosses the reth 
level through Xn (The root need not be labelled). Hence, the set dr of ends of r is naturally 
embedded in the sequence space B^ . 



2.4 Embedded Labelled Galton-Watson Trees Tr 

The construction of embedded Galton-Watson trees in the anisotropic contact process 
resembles the ones for anisotropic branching random walks and isotropic contact pro- 
cesses [15|. The offspring of a vertex x £ Tr will be vertices y at distance r from x such 
that there is a downward infection trail beginning at x and ending at y. It is apparent that 
the concatenation of infinitely many such downward infection trails gives rise to an end of 
Tr that is contained in the limit set A. 

Definition of r^. Fix an integer r > 1. Define generations Vn{r) = Vn G C^r inductively 
as follows: 

(i) Vo = {1}. 

(ii) For each x £ Vn, the offspring of x are those y S J~-nr+r such that there is a downward 
infection trail from x to y, beginning at time Sn, the first time when x is reached, 
and first reaching y at time 5^+1 (suppose that 5o = 0), and this infection trail is 
the first downward infection trail beginning at {Sn, x) to reach y. The offspring vertex 
y has parent x. 

(iii) Vn+i is the set of all offspring of vertices in Vn- 

Note that the random times Si are stopping times. For each vertex y 7^ 1 of the tree r^, 
define the label of y to be the word x'^y, where x is the parent of y. Each label is a reduced 
word of length r which ends in the letter a. Thus, the set of labels is finite. 

Lemma 1 x,. is a labelled Galton-Watson tree. 
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Proof. The infection trails associated with different x € Vn are distinct and the existence 
of an infection trail depends only on the percolation structure in T(x) after time 5„, 
thus, is independent of the pre-Cnr history of the contact process initiated with Hence, 
the trails S,x, x S Vn, involve nonoverlapping regions of the percolation structure and do not 
overlap the region of the percolation structure that determines Vn, whence by the strong 
Markov property, are mutually independent. Since each vertex of Vn has final letter a, the 
contact processes initiated by the trails ^j,, x G V^, are all "oriented" the same way, relative 
to the level structure of the tree T* . This implies that the offspring distributions for vertices 
X £ Vn as probability distributions on the set of labels are all the same as that one of the 
initial vertex 1 in Vq. As a consequence, is a labelled Galton- Watson tree. □ 

Observe that the relation between the metrics d'^ for the tree r,. and for the full tree 
r is 

d^j{x,y) = dar{x,y). (2.3) 
The mean offspring number //^ for the Galton- Watson trees is, by construction, 

x€C*. 

2.5 Geometric Decay and Strict Monotonicity of 

It is not difficult to show that, in the isotropic case, A G TZ^, implies /3(A) = 1. The reverse 
direction of this statement is more subtle but has been shown in ||l^, more precisely, if 
A ^ TZs, then /3(A) < l/y/2d — 1, with strict inequality holding for each A in the interior 
of the weak survival phase and equality emerging at the transition. The feature /? < 1 
guarantees exponential decay of the infection probabilities u„ in n, for every integer n > 1. 
It is essential to verifying both, the strict monotonicity and continuity of the relevant growth 
variables of the contact process in the infection parameter. 

In fact, for the anisotropic contact process as well, /3i(A) < 1 for each i £ A and for 
A 7^3 would be a desirable property in proving strict monotonicity and continuity of 
certain functions in the infection parameters. Different values of these functions at TZi DTZI 
and TZ^ n TZ^ could then be used to argue that the region 7^2 has nonempty interior. To 
help the understanding of the rest of the paper, we here outline our strategy in establishing 
the existence of a weak survival region. We shall prove the following: 

Strategy of proof. 

(A) If each /?j(A) < 1, then we have continuity and strict monotonicity of certain functions, 

(B) a set Pc of discontinuities is located, 

(C) as A with each Pi{X) < 1 approaches any point in Dc, each /3j stays bounded away 
from 1. 

Prom (A) — (C) together with several other considerations, it will then follow that in 
fact the set Vc coincides with the boundary TZ^ n TZs and that (TZ^ n 7?.i) n {TZ^, n 7^3 ) = 0. 
It will follow that Vc separates the regions TZi U 7^2 and 7^3, and that the functions at 
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hand are strictly increasing and continuous for A G int(7^i U 7^2) and left-continuous for 
X G TZi U7^2 such that each > for d > 1. Sections 5.5 and 5.6 are the principal sections 
in our discussion of identifying the three regions and classifying their boundaries, whereas 
the sections beforehand introduce machinery and give a full account on the specifics about 
the functions which are to be continuous and strictly monotone. Some of the latter parts 
may appear a bit technical. 

The distinguished set of all A that have each /?i(A) < 1 is given a name as follows. If 
Q = {A G R"^ : each Afc > 0}, define the sets 

JC = {A G Q : /3i(A) < 1 for each i £ A}, 
/Co = {A G /C : Aj > for each i G A}. 

Geometric decay of Ux- We begin with some estimates on the decay of the infection 
probabilities. We point out that the upper bound for u^(A) in ( |2.5| ) below is uniform in A 
for d > 2, that is, if we assume that there is some e > such that at least two infection 
rates A,, \j > e for two distinct i,j G A-\-. If A is restricted to those A G /C, then it can be 
shown that there is some p = p(e) < 1 such that each Pi{X) < p < l- 

Lemma 2 For 6V6Ty fixed A, each ijitegeT ^ 0, dfid evevy x G Qn: 

u,{X) < [max/5,(A)]^ (2.4) 

aeA 

For each A G /C and each x € Qn, there is some constant < 7 < 1 such that 

n,.(A) < 7". (2.5) 



Proof. The proof proceeds by induction over the distance of the vertices from the root. Fix 
A. Without loss of generality, we may assume that A;, > Aj for every i £ A. Let Un = bb . . . b 
with lunl = n. For n = 1, clearly, Uyi(A) = n;,(A) > Ui{X) for every i £ Ahy the definition of 
the functions n^(A) and the one of the rules of infection of the contact process. Now assume 
that Uy^^ (A) > Ux{X) for every x G Qk and every A; = 1, 2, . . . , n — 1. Thus, no x G Gn-i has 
larger probability than y„_i to ever be infected. But upon infection of Un-i, no vertex in Gn 
has larger probability to be infected than y„ because of the rules of infection of the contact 
process, the homogeneity of the process at each vertex, and the fact that, along any path 
bbb...., the neighbourhoods of the vertices look the same. Hence, Uy^{X) > Ux{X) for every 
X £ Qn- Since this argument is valid for every integer n > 0, the desired results follow from 
the fact that Uy^{X) < /3{,(A)" for every n > 0. □ 

Strict Monotonicity Properties. For each infection parameter A G int(/C), there 
exist directions of increase and decrease, respectively, where the infection probabilities Ux 
strictly increase or decrease, respectively. It is not apparent, however, whether there is an 
easy criterion to decide for each vector pointing away from some A, whether the infection 
probabilities have strict monotonicity properties. We will construct directions of strict 
monotonicity for each A by "thinning" the percolation structure of the contact process. This 
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useful idea to modify the percolation structure originates from [13|, where in the isotropic 
case strict monotonicity of the function f3{-) in the infection parameter was shown. In the 
anisotropic case, the precise statement is more subtle, though. For this purpose, we first 
describe a modified contact process that arises when the percolation structure is manipulated 
by a set of Bernoulli-p random variables, where p is precisely chosen to decrease the smallest 
infection rate Ac, say, to \c < \c- The next proof will show that tuning the smallest infection 
rate, forces all infection rates to decrease, which might possibly be more than necessary but 
it will suffice to lead to strict monotonicity of many characteristics of the contact process 
that will be discussed later. 

Suppose that Ac = min^e^ Xa is the smallest positive infection rate of the contact process 
At with infection parameter A and recovery rate 1. Fix Ac < Ac- To each infection arrow co 
of the percolation structure, there is attached a Bernoulli-p random variable ^(^. These are 
conditionally independent, given the realization of the percolation structure. Choose the 
value p = P{Cuj = 1} so that 

Ac = pAc/(l + qXc)- 

where g = l— pE(0,l). Call this assignment of a set of Bernoulli-p random variables to the 
percolation structure a. p-thinning of the percolation structure relative to Ac- A version A[ of 
the contact process with infection rates pXi and with recovery rates rj = l + gAj, i G A, may 
be constructed using the augmented percolation structure by (1) first setting up a modified 
percolation structure by changing every infection arrow lv such that = to a recovery 
mark *, then (2) defining A'^ to be the set of all vertices y for which there exists a directed 
path from (1,0) to {y,t) in the above manipulated percolation structure that does not pass 
through any recovery marks *. The new contact process A'^ is a time-changed version of a 
contact process with infection rates A' = i{X'a}aGA+) recovery rates r' = ({?'a}ae^+) 
with A^ = pXi/{l + qXc) < Aj and = (1 + qXi)/{l + gAc). Thus, since r'^ = 1 and > 1 
for each i £ A, there is a contact process A^ that has the same percolation structure as A'^, 
whose recovery rate is 1 and whose infection rates are no larger than A^. 

Note that it is enough to state the next result for a contact process that has only positive 
infection rates Aa (because otherwise the set of generators A+ may be updated). 

Proposition 1 Let At be a contact process with infection parameter X and recovery rate 1. 
Let A G /Co with < Ac = min^g^Aa. Assume that its percolation structure is p-thinned 
relative to Ac for some < p < 1 so that a contact process At is obtained with infection 
parameter X and recovery rate 1 with < Ac = pXc/{l + qXc) < Ac and Aj < pXi/{l + qXc) < 
Xi for every i £ A. Then there is some constant < uj < 1 such that for every n and x £ Gn, 

u,{X) < n,(A)u;". (2.6) 

In particular, (3x{X) < Px{X) for every x £ G, where (3x is defined in ( ji.^ j. 



Proof. The ideas of the proof are essentially the same as those given in |13|, Proposition 9, 
for the isotropic case. We outline the necessary modifications and refrain from repeating the 
entire (not short) proof. The argument rests on the previously described modification of the 
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percolation structure of the contact process by p-thinning by means of Bernoulli-p random 
variables relative to Ac- Choose the unique value p so that Ac = pXc/ (l+^Ac) . Assume that At 
is a contact process with infection parameter A, and recovery rate 1 and that is a contact 
process with infection parameter A and recovery rate 1 with Ac = A^ = pAc/(l + qXc) < Ac 
and Aj < A- = pXi/{l + qXc) < Aj for every i £ A. Thus, all infection rates have been 
decreased by the thinning process. 

To verify (^.61) , it suffices to show that there is some uj G (0, 1) such that the contact 
process A'f-, gotten from the p-thinned percolation structure, with infection parameter A' 
and recovery rates r' has Ux{X') < Ux{X)uj^ for every x £ Qn and every integer n > 0, 
because the contact process A^ is equivalent to A[ in the sense that they have the same 
percolation structure, and thus, Ux[X) = Ux{X') for every x £ Q. 

Let Gx be the event that the contact process At infects vertex x € Qn a^t some finite time, 
and let G'x be the corresponding event for the contact process A[, thus, P{Gx} = Ux{X) 
and P{G'x} = Ux{X') (Here we abuse the notation since except for this proof, Ux refers to 
a contact process with recovery rate 1). By construction, G'^ C Gx because every infection 
arrow in the modified percolation structure occurs in the unmodified one as well, and on 
the other hand, every recovery mark * in the unmodified percolation structure is retained 
in the modified percolation structure. On the event Gx, there is at least one directed path 
through the unmodified percolation structure that leads from (1,0) to x x (0, oo), in fact, 
there may be many overlapping such paths. Call an infection arrow to in the unmodified 
percolation structure essential for the event Gx if (1) event Gx occurs, and (2) changing 
to from an infection arrow to a recovery mark would destroy all directed paths from (1,0) 
to X X (0,oo). Define A^^^ to be the number of essential arrows for the event Gx, when Gx 
happens, and Nx = when Gx does not happen. Since removing any one of the essential 
arrows would disconnect (1,0) from x x (0, oo), in order that event G'^ occur it is necessary 
that Nx > 1 and that = 1 for every essential arrow uj. This event has conditional 
probability p^'' , given a realization of the unmodified percolation structure. Hence, 

Ux{X') = P{G'J = Ep''^lG^, 

where /|.} denotes the indicator function. If it were the case that Nx > cn on Gx for 
some positive constant c, then it would follow that Ux{X') < Ux{X)p'^'^. On the other hand, 
if for some p < 1 \i were the case that P{Nx < cn\Gx} < p"", then it would follow that 
UxiX') < Ux{X){p^'^ + p"). Therefore in both cases, it would follow that there is some 
< UJ < 1 such that u^(A') < Ux{X)uj'^ and, since this argument holds for every n > and 
X G Gn, prove ( |2.6D and, by subadditivity, that Px{X) < /3x{X). Now, verifying Lemma 3 in 
|l3| (there stated in the isotropic case) for the anisotropic contact process will complete the 
proof: 

There exist constants c < oo and < p < 1 such that for all sufficiently large n, 

P{Nx < cn\Gx} < p" 

( p3|| , Lemma 3). Note that this statement refers only to the unmodified percolation struc- 



ture. In the proof of Lemma 3 given in [13|, only features common to both the anisotropic 
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and isotropic contact processes are relied on except for a single passage in the very last 
paragraph of the paper, which requires a concrete estimate of the underlying anisotropic 
contact process, that is, an upper bound 7", < 7 < 1, for the probability n^(A) for every 
sufficiently large n and each x ^ Qn- But, since A G /Cq C /C, this result follows from ( p.5D 
in Lemma This finishes our proof. □ 

We say that the directions given by the vectors A — A and A — A are directions of increase 
and decrease for A, respectively, for the contact process with infection parameter A (and 
recovery rate 1), where A is defined in Proposition |l] for any < p < 1. Note that any 
function in A that is strictly increasing along directions of increase for A necessarily is 
strictly increasing if every component of A is increased. 



3 Potential Functions and 77 < 1 

The proofs of the principal results of this paper rely on classical results from the theory 
on Gibbs states as developed in ||3|, |2^. They enable us to "normalize" the ensemble of 
infection probabilities Ux so that, for any e > and at all fixed sufficiently large distances n 
from the root vertex, there exists a shift-invariant probability measure concentrated on the 



set of vertices in with n(l — e) < k < n(l + e) (see e.g. the remarks surrounding (3.24) 
through l^?rj\)). 



We begin with describing this approach via potential functions, quote a result that 
applies the theory to some matrices that satisfy a Holder condition, and proceed to apply 
these tools to the collection of infection probabilities for the contact process. Our goal is 
to define Gibbs states supported by the geometric boundary fi, which will allow us to keep 
track of the growth of the number of infection trails that wander off to infinity and of those 
that return to the root vertex at some finite time. From the properties of Gibbs states, we 
will derive estimates that describe the dispersal of the contact process on T in space-time, 
as for instance, the rightmost and leftmost positions of the infection at time t, which both 
will turn out, almost surely on the event of survival, to move at a distance from the root 
that is linear in time. 



3.1 Background: Thermodynamic Formalism 

Recall that A denotes the set of generators of Q and their inverses. Let = {one-sided 
infinite sequences from and A!^ = {two-sided infinite sequences from A^. Recall S to 
be the set of all doubly infinite reduced words ^ = {xn)'^=-ao from A (reduced means 
that, for every n, Xn+i 7^ The subset S is a closed subset of the coding space A^, 

thus, is compact in the metric d^. Let a : S — > S denote the forward shift on S, that 
is, a{xkXk+i ■ ■ ■) = Xk+iXk+2 ■ ■ ■ for every k. Note that a is Lipschitz continuous. For any 
function / : S — > R, define 

Snf = f + foa + foa^ + ... + fo 

for every n > 1. 
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Potential Functions, Gibbs States and Pressure. Define the n-cylinder sets 



r(i) = r(iii2 ...in)={^ = (xfc)^.^ G S : Xj = ij, 1 < j < n} 

for every n and i = iii2 ■ ■ ■ in ^ Gn- spirit of Q, for any Holder continuous function / on 
S, there is a real constant P{f), some constants < Ci < C2 < cxd, and a unique a-invariant 
probability measure /ij on the Borel sets of S such that for each i — i\i2 ■ ■ ■ in ^ Gm 

- exp{5./(,) - nP(/)} ^ ^^-^^ 

for every j G r(i). The measure /i/ is called the Giftfts stoie with potential function /, and 
the constant P{f) is called the pressure of /. 

We list some features of the pressure function that will be useful to our subsequent 
analysis. Two Holder continuous functions / and g are called cohomologous if there exists a 
Holder continuous function h such that f — g = h — hoa. If two functions are cohomologous, 
then they have the same pressure and the same Gibbs state. The function a Pio-f) is 
continuous ||2^. For every Holder continuous / and for every integer n > 1, the pressure 
functional satisfies P{Snf) = nP{f). Most importantly, if / < 0, then there exists a unique 
constant 5 > such that 

P{5f) = 0. (3.2) 

The pressure and the measure-theoretic entropy of the Gibbs state are related to each other 
by the Variational Principle (Q, Theorem 1.22) since prominently the Gibbs state for tf, 
any t > 0, is the unique equilibrium state for —tf. Choosing t to be the 5 that nullifies the 
pressure yields 

h{^^sf) = P{5f) - j 5f dfisf = -J^f dm- (3-3) 

Counting problems. For < C < Ij define 
00 

■A*{C) =\J{i^Gn: Snfii) < log C and Suf{i) > log C, < n} (3.4) 

n=l 

to be the set of comparable finite sequences relative to / at scale C- The set A*{C) consists 
of those finite sequences of possibly different lengths n such that Snf takes a value just 
below logC. Observe that, since / is bounded, for any i S A*{C)j Snfii) differs from logC 
by at most ||/||oo < cxd. Note that for every sequence i G S, there exists a unique n such 
that the finite sequence iii2---in is an element of A*{C)- The following result, borrowed 
from , is presented along with a proof because the latter offers some insights that will be 
of use in the subsequent discussion. 

Proposition 2 (Hueter and Lalley Q, Proposition 2.1) Let 6 > be the unique pos- 
itive number such that P{Sf) = 0. Then there are some suitable positive finite constants C3 
and C4 so that, as C, ^ 0, the cardinality of A*{C) i^ given by 

CsC' < \A*iC)\ < C^C'- (3.5) 
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Proof. Each element i G ^*(C) gives rise to a cylinder set T{i). Moreover, the cylinder 
sets {T{i) : i G A*{C)} are pairwise disjoint and their union is the entire sequence space S. 
Consequently, 

By the definition of 5 and the Gibbs state, there are constants < Ci < C2 < oo such that, 
for every i G A*{C), the measure of the cylinder set T{i) satisfies 

Ciexp{6Snf{i)} < ^5/(r(i)) < C2exp{6Snf{i)}, 

where n is the length of i. But by the defining property of A*{C), there is a positive constant 
C3 < 1 such that for every i £ A*{C)j 

C3C < exp{5„/(i)} < C- 

Assembling the last three displayed formulae yields 

C4 E c'<i<c2 E 

for a suitable constant C4, which proves the advertized inequalities. □ 

"Most" sequences in A*{C) approximately /^^/-distributed. Indeed, by the Birkhoff 
ergodic theorem, for every Holder continuous function (7 : E — > R and every e > 0, 

\{i = hi2...in£A*{C) ■■ max \^ME^ _ t [ gdfisf\>e}\=o{C^). 

0<t<i n J 

Therefore, if we let = logC/ / gd^sf, then "most" sequences in ^*(C) have lengths be- 
tween n^(l — e) and n^(l + e). These observations indicate that the set A*{C) is nearly a 
set of sequences of length that are approximately "generic" for the measure fisf- In light 
of the Shannon-McMillan-Breiman theorem, the cardinality of the latter is approximately 
exp{/i^^^n^}, where /i^^^ denotes the entropy of the measure fisf, whereas in view of Propo- 
sition the cardinality of the former set is of the order (^~^. This confirms the variational 
principle in ( |3.3D , which implies that 



Next, we state a result of |11], Proposition 5.1, an extension of the Perron- Probenius 
theorem. For any matrix M, let ||M|| = sup„^Q(|Mt;|/|f |) denote the usual matrix norm. 



Theorem 9 (Lalley |11|, Proposition 5.1) Assume that Mx is a nonnegative aperiodic 
nxn matrix andx is a Holder continuous function (with some exponent) onA'^. Then 

there exist constants C < 00 and < a < 1 and Holder continuous functions 7 : S ^ R 
and v,w : T, ^ V+ = S R" : X]?=i Ui = 1 and Ui > for every 1 < i < n} such that for 
every ^ = {xk)^^_^ G S and every integer n > 0, 

||e-^"'^«)M,,M,., . . . - i{a'-i)v{i)w{a''i)'\\ < Ca^, (3.6) 
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where 



7(0 = l/w{0'v{0, (3.7) 
M.,M<^0 = e'^(«)r;(0, (3.8) 



and 



w{a-^0'M^, = e^(-0 ^(^)*. (3.9) 

Both v{(^) and ip{S,) are functions of the 'forward" coordinates xi,X2,--- and w{£,) is a 
function only of the "backward" coordinates . . . ,X-i,xq. 

Thus, the Perron-Frobenius theorem has x a constant function. We will consider 

the special case when the matrices Mj, are 2d x 2d, thus, the size of the matrices being 
fixed and independent of the number of terms in the matrix product in (|3.6|). The proof 



in |11| carries over unmodified. Observe that, for each n, v{^) and w{^) are right and left 
eigenvectors, respectively, of M^^Mrj . . . M^^ with associated eigenvalue exp{Sn<f{0} ■ 



3.2 Critical Exponent 

Recall that Gm denotes the set of all vertices at distance m from the root vertex. For every 
A, define 

00 

Tu = ruiX) = inf{r > : ^ ^ u^iXy < 00}. (3.10) 

m=0 xeGm 

This exponent r„(A) takes some finite value for A € int(/C) and is nondecreasing in each Aj 
because Ux{X) is nondecreasing in each Xj. For instance, it is easy to see that, for A T^i, 
we must have ^^(A) > 1. Some intuitive values yet subtle to prove are as follows: we will 
ultimately show (Corollary |ll|) that, for d > 1, we obtain ru{X) = 1 for A G TZi Ci IZ2 and 
fu{X) = 2 for A S 7^2 n 7^3. We begin to view some easy but crucial facts. 

Lemma 3 For each t > ru{X), we have 

00 

m=0 X&Gm 



Proof. This is readily concluded from the inequality < u^". □ 

Proposition 3 Let A G /Co with < Ac = miiiaeA Xa cmd let X he an infection parameter 
that corresponds to the p-thinned percolation structure for some < p < 1 relative to Ac 
such that < Ac = pAc/(l + qXc) < Ac. Thus, X lies in a direction of decrease for X. Then 

ruCX) < ru{X). 
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Proof. Fix some < p < 1, let A € /Co, and let A be a parameter that corresponds to the 
infection parameter gotten from the p-thinned percolation structure relative to Ac in the 
sense that the percolation structure is preserved. Choose A such that the contact process 
has recovery rate 1. In view of Proposition there is some constant < a; < 1 such that 
for every n and x G Gn, u^^X) < Ux{X)io"'. Moreover by Lemma ^, for each t > r„(A), we 
have Em=oExeg,„4(^) < ^- Consider 



oo 



^ ^ u,(ArW-^ < ^ ^ (u;™rw-^n,.(Arw-^ 

oo 

= E E [w'"^""^^^"'V^*x(A)'] u,(A)''"W-^+^ (3.11) 

m=0 x&Qm 

for some reals s,5 > 0. In view of parallel arguments, resting on the minimal positive 
infection rate Ac, to the ones employed in the proof of Lemma ^ for the maximal infection 
rate along with ( |1.5D , for any e > 0, there exists some such that for every m > N and 
X e Gm, we have > (1 - e)-[min,e^ = (1 - e)'"(/5c)'" > (/3c/2)'" with > 
because Ac > 0. Since w < 1, for every s < r„(A), we can choose 6 small enough so that 

for some < 7 < 1 and sufficiently large m. In addition, for each s < 5, we obtain 
J2m=QJ2xegm'^^^^y'^^'^''^'^^^ ^ ^ some positive finite constant C. Observe that there 
exist sufficiently small s,5 > that can satisfy the former and latter conditions. Hence, 
all sufficiently advanced terms in the summation X]m=o righthand side of ( 3.11| ) are 



bounded above by 7"*C for some positive finite constant C. Therefore, summation over a 
geometric series leads us to conclude that there is some positive s such that 



00 



^ J2 n..(A)'^"(^)-^ < 00. 
m=o xeg,n 

Consequently, r„(A) < r^iX), as desired. □ 

Proposition 4 Suppose that A' G /C Pi /C'^ and that A lies in the interior int{lC) and in a 
direction of decrease for A'. Then 

ru{X) < r,(A'). (3.12) 

Proof. Assume that A lies in the interior int(/C). Then there exists a ball B in R'^, centered 
at A which is completely inside /C. Since each ball contains a suitable multiple of each 
unit vector in R'^, for every < p < 1, the ball B contains a line segment Lp in the 
direction of decrease for A, equivalently, increase that contains A, defined by the p-thinning 
of the percolation structure relative to the minimal infection rate. Hence by Proposition 
each point A^, of Lp that has strictly smaller distance to the boundary IC D IC'^ than A has 
TuW < ru(A*). But ruiK) < r„(A') for A' G ^ fl Hence, r„(A) < ru(A*) < r„(A'), as 
advertised. □ 
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3.3 The Infection Probabilities satisfy a Holder Condition 



Perhaps the most interesting aspect of this and the next section is how to prove a Holder 
condition for the cohection of infection probabihties in /Co so that the existence of a Gibbs 
state is guaranteed for the contact process and the theory described in Section 3.1 comes 
to fruit. 

Proposition 5 There are constants < 7 < 1 and < C < 00 such that for every 
A G int{ICo), every p > ru(A), every integer k > 0, and for every sufficiently large integer n, 



Proof. Fix A € int(/Co)- Assume that Af, > Aj for every i £ A. Let denote the probabihty 
for the isotropic contact process with infection parameter e > and initial infection at the 
root vertex that between time and 1 there is no recovery mark * at the root 1 and that, 
for some i £ A, the root vertex infects its neighbour i. Thus, = e~^(l — e"*^) > 0. Fix e 
such that each Aj > e. Then for every x G ^fc, we obtain 

Fix some integers k,M > 0. Since p > ru(A) = r„, we can choose some 5 > such that 
p — r^ — 5 > 0. As a consequence, by Lemma ||, for each x £ every n > kM, and each 
y e Gn, 

Since, by choosing M sufficiently large, (3i,{\Y~'^'^~^ / {a^Y^^ < 7 < 1 for some 7, and thus, 
the righthand side of the last display is bounded above by 7". Combining this with Lemma 
^ yields, for every kM < n, 

E^ 'UxiX2...X„{X) _ {UxiX2...Xn{X))^ ,,-.r^+S 

< i^c 

for some positive finite constant C. This ends the proof. □ 



3.4 Transition Matrices and Potential Functions 

The Holder continuity allows us to define a Gibbs state by means of certain matrices. 
Whereas there are numerous reasonable choices of definition for transition matrices in the 
context of the contact process on T, the one we propose is natural and has some nice 
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interpretations for some exponents. If nothing else is said, we shall assume that each 
Afc > and that the generators with zero infection rates have been eliminated from A-\-. 
For each xiX2 ■ ■ ■ Xk-i G Gk-i and for all i,jGA with i ^ 2;^:*;^, define 

Tij = T{xi . ..Xk-l,i,j) = {Xk+l . ..Xn-l G Gn~k-1 ■ Xk+1 / Xn~l / j'^} (3.14) 

and let !F{xi . . . Xk-i,i,j) be the empty set if i = x'^^^-^^. Then, for each real p > 0, 
all integers n — 1 > k > 1, and each xiX2 ■ ■ ■ x^-i S Gk-it define the 2d x 2d matrix 
Hp{n;xiX2...Xk-i;X) by 

iH,in;x^X2...Xk-r,X)h= ^ ( ^-^-^•■-"(^) V (3.15) 



Xk+lXk+2 ■ ■ ■ Xn-l G J^ij 



for all i,j G A, thus, equal zero for i = x'f^\. Note that for p = 1, k = 1, the (z,j)-entry of 
Hi{n; 1; A) equals the expected number of vertices at distance n from the root 1 that are 
ever to be infected whose word representation begins in i and ends in j. Similarly, for p = 2, 
k = 1, the (z, j)-entry of H2{n; 1; A) is a lower bound for the expected number of vertices 
at distance n from the root whose word representation begins in i and ends in j that are 
ever infected and upon infection are ever to send the infection back to the root. Moreover 
for each real p > 0, l^Hp{n; 1; A)l = Z]xge„ where 1 denotes the 2d-vector of all ones. 
Importantly, we recover 

l*/7i(n;l;A)l = ^ u^, (3.16) 
l'H2{n-l;\)l = <■ 

X&Qn 

By subadditivity, it is an easy observation that 

\Hp{n + k-l; xiX2 . . . Xk-i, A)|| > \Hp{n] 1; A)|| 

for every /c > 1, in other words, there might be positive probability that, by the time of 
infection of vertex xi . . -Xk-i, some other infected vertices line up in the complement of 
T{xi . . -Xk-i) to reinforce the infection. Observe that for every A G /Co and each p > 0, 
the matrix Hp{n;x;X) is an aperiodic, irreducible, and nonnegative matrix, thus, a Perron- 
Frobenius matrix. Therefore, the Perron-Frobenius theorem tells us that Hp{n; x; X) has 
a largest positive eigenvalue. We will see in Proposition ^ below that the lead eigen- 
value equation can be solved explicitly. In light of Proposition ^, for p > ru{X), the 
matrix Hp{n;x;-) satisfies a Holder condition. The mapping x — > Hp{n;x;-) is given by 
X1X2 ■ . . Xk-i Hp{n; X1X2 . ■ ■ Xk-i, •)• 

Proposition 6 For every X G ini(/Co)) CLnd each real p > ru{X), the matrix x Hpin; x; •) 
is Holder continuous (with some exponent). 



Proof. Fix such A and any p > ru{X). Since the matrix norm ||j4 — of the difference of 
two nonnegative matrices A and B is bounded above by the maximum max(||j4||, ||-B||) of the 
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individual matrix norms, and in turn, the norm of any matrix is bounded by the maxi- 
mum of the absolute values of the ^-matrix entries times its dimension, it suffices to verify 
that the maximal absolute values of the entries of the differences Hp{n; xiX2 ■ ■ ■ Xk-i] •) — 
Hp{n; X1X2 . . . Xk+i-i] ■) and Hp{n; X\X2 ■ ■ ■ Xk-i; •) — Hp{n + I; X1X2 ■ ■ ■ Xk-V, •) are bounded 
above by C'j"' for some constants < 7 < 1 and < C < 00, for all integers k,l > and 
for every sufficiently large n. But it is an elementary exercise by the results in Proposition 
P to show that each entry of the latter four matrices is bounded above by C7", where C 
and 7 are independent of k,l, and n. This finishes our proof. □ 

Recall that S denotes the set of doubly infinite reduced words from A. Notice that 
in the subsequent discussion, we shall not always be careful to distinguish and S. Now 
we are ready to define the potential function ipp{x) = ipp-\{x) for every p > r„(A) by 
applying Theorem |9| with M^^ ■ ■ ■ Mx^_^ = Hp{n; xi . . . x^-i; ■), V = v and W = w with 
V,W :T. = {u e K'^'^ : Y?iii Ui = 1 and Ui > for every l<i<2d}. Then for every 

X = . . . X1X2 ... G S, every integer k > 0, and sufficiently large n, define 

Hp{n; XIX2 ■ ■ ■ Xk-i; A)y(cj"-'=x) = e^^-k'Ppi'')v{x). (3.17) 

A similar equation represents the left eigenvector W{x) associated with Hp(n; xi . . . x^-i; A) 
(for more details, see Proposition 5.1, [|ll|). The Holder continuity of (fp = ipp-x immediately 
follows from the Holder continuity of V and the map x — > Hp(n; xi . . . Xk-i] A). To facilitate 
notation, we will often just write ipp. Observe that H p[n; X1X2 ■ ■ ■ Xk-i] X) has eigenvalue 
exp{5„_fcV3p(x)} with corresponding right and left eigenvectors V{x) and W{x). Also, note 
that for p > r„(A), we have Lpp < 0. Combining ( |3.6| ) with definitions ( |3.15[ ) and (3.17) 
yields 



r, ,-, T, , PG \ ^2:1X2. ..Xfc (A) 



Cexp{SnM^)}{l + 0{a^)) (3.18) 



for some < a < 1, for every x = . . . X1X2 ... G S, all integers k and n, and every p > ru(A), 
where the implicit bound in the 0{a'') term is uniform in x. The constant C may be bounded 
by < Ci < C < C2 < cxD with Cj independent of x, k and n. 

For every x G S, the functions Sn^p{x), and thus, Hp{n + k;xi. . -Xk-i; ■), depend on 
the relative frequencies of the generators in the reduced word X1X2 ■ ■ -Xk-i and the order 
of the letters as well. 

Matrix Entries 6j. Next we regard to express each matrix H2{n]XiX2 ■ ■ -Xk-i^X) as 
the (n — A;)-th power of a Perron-Frobenius matrix -B2,n,xi...a;fc_i (A) which (a) preserves the 
lead eigenvalue, and (b) has a certain form that specifies the allowed transitions as the 
matrix does in the following result. 

Proposition 7 For any positive number p, define the matrix Ap to be the 2d x 2d matrix, 
indexed by elements of A, whose entries are given by 

{Ap)ij = a'. ifj^i-\ (3.19) 
= ifj = i-^ 
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for some real numbers aj > 0. Assume that Uj = aj^i for each j £ A and that Ap is a 
Perron- Frobenius matrix. Then the lead eigenvalue of Ap is the unique positive solution ap 
of the equation 



a' 



E ^ = 1- (3.20) 

Moreover, for each p > 0, the function strictly increasing for each i £ A. 



Proof. Identity ( ^.20 ) is shown by rearranging the eigenvalue equation ApU = UpU. If 
Up is the lead eigenvalue, then the vector u must have nonnegative entries, not all zero. 
Let < s = X^ieyl ■ This identity combined with the eigenvalue equation UpUj = 
Y.i<aA ^i'^i ~ ^j'^j yields 

Uj = s/{ap + ap. 

Multiplying both sides by adding over all i £ A, and dividing both sides by s gives 



( 3.20| ). In addition, the strict monotonicity of Op in each aj is an immediate consequence of 



the form of the equation ( 3.20| ) and the facts that the map x x/{ap + x) is continuously 



differentiable and has strictly positive derivative for positive x. □ 

Now, for any positive number p and A G /C, define a2dx2d matrix Bp = Bpn,xi...Xk-i 
indexed by elements of A, so that (a) their entries bj{X)P = bj(X, xiX2 ■ ■ ■ x^-i)^ > come 
in the form ( ^.191) , that is, 

{Bp{X)),, = b,{Xr if i/r\ (3.21) 

= if j = 

and (b) the bj = bj{X) satisfy ( ^.20D with 

= b"^,, p = 2, and a2 = exp{Sn-W2{x)/{n - k)}. (3.22) 



Hence, the bj, j S A, link (f2 and the matrices Bp. In this construction, the power p in (3.22) 
needs to be chosen sufficiently large to guarantee the existence of fp and to assure that the 
bj are well defined. This is possible since one shows that r„(A) < oo for A G /C. Our choice 
p = 2 is the smallest possible because it will turn out later (Corollary ^) that ru(X) < 2 for 
A G int(/C). Thus, i?2~^ has lead eigenvalue exp{Sn-k'^2ix)} . Both Bi and B2 are Perron- 
Frobenius matrices. The entries bj of Bi and bj of B2, respectively, depend on n and the 
sequence X1X2 . . .x^-i. We shall not always be careful to indicate these dependencies and 
assume that they are clear from the context. Observe that, if /3i(A) = 1, then for each 
X € Gk-ij we can choose bi{X,x) = 1, due to the construction of the matrix Bi, equation 



( 3.20| ), and the fact that Uy is independent of the choice of the beginning finite segment 



X1X2 . . . Xk-i of an infinite word y = . . . x\X2 ■ ■ ■ ■ An analogous argument shows that, if 
/3i(A) = 0, we may choose 6j(A, x) = for each x £ Q. Oi course, the converse also holds. 
The following observation is a summary. 

Lemma 4 Ifbi{X,x) < 1 for each x £ Q and i £ A, then X £ fC. Moreover, bi{X,x) > for 
each X £ Q and i £ A if and only if f3i{X) > 0. 
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The next result shows that the functional Sn~k'-P2 is compatible with the construction 
of the matrix i?2- 

Proposition 8 For every sufficiently large n — k > 0, the functional Sn-k^2 is additive, 
that is, the function (/?2 : 5^ — > R satisfies 



where a denotes the forward shift on S. 



Proof. Let x = . . . xiX2 . . . G S. Since 02^^ is the lead eigenvalue of H2{n; x\X2 ■ ■ ■ a^fc-i; A) 
and i?2~^' thus, the right eigenvector u of B^~^ is dominated by a positive finite 
constant multiple of the right eigenvector of Ji2{n\x\X2 ■ ■ - x^^^i^X) and vice versa, the 
(right) eigenvectors of B^'^ and B.2(n\ x\X2 ■ ■ ■ Xk~i; A) differ by at most a positive finite 
constant multiple. Additionally, in the eigenvalue equation B2~^u = a2~^u, the same 
eigenvector u is relied on for every n — k. Hence, iteration brings 

exp{Sn-kV'2{x)}V{x)ci = H2{n;xiX2...Xk-i;X)V{a"-~''x)ci 

< B'^'^u 

= a2B^-~^-^u 

— ^,2 Dn-fc-2,, 

= a2 -D2 ^ 
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n~k—l 

= exp{ ^ ip2{cF''x)}u. 

Similarly, the reverse inequality holds if the constant ci is replaced by a constant C2- Since 
for sufficiently large n, the multipliers Cj are negligible as compared to the argument of the 
exponential function, this completes the proof. □ 



3.5 Invariant Measures 

Next we recollect some facts from Section 3.1. An elementary exercise shows that r„(A) < 00 
for A G /C. We have seen that ipp < for each p > r„(A). Thus, there is a unique shift- 
invariant probability measure fi^p^ on the Borel sets of S such that there is some real constant 
P{(Pp) and some constants < C3 < C4 < 00 such that, for each m, i G Q^, and cylinder 
set r(z), we have 

C3exp{5m<^p(j) - mP{Lpp)} < p^p{T{i)) < C4 exp{S'm93p(j) - mP{Lpp)} (3.23) 

for every j G T{i). From (p.4|), recall 

00 

^p(C) =\J{i^Qn: SnVp{i) < logC and SkiPp(i) > logC, VA; < n} (3.24) 

n=l 
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for every 1 > C > and that, for p > r„(A), the sum J2ieA''{0 t^^'PpC^ (''')) ~ ^^^^^ 
6 = 6p > is such that the pressure P{6ipp) = 0. Define the expectation of ipp with respect 
to the probability measure jj^Sipp by 



'Pp = J "Pp^f^Svp- (3-25) 
The Birkhoff ergodic theorem imphes that, for every e > 0, 

\{i = hi2 . . . i„ e A;{0 : max _ tTp^\ > ,}| = oiC'), (3.26) 

thus, "most" sequences in A*p{C) have lengths between n^(l — e) and n^(l + e), where 

nc = (3.27) 

Next, consider the limit points uj = ujiuj2 ... € in the weak survival phase. Define the 
distribution /i„ under P of the process a;„,C(;„+i, . . . , that is, for any Borel set V C ^l, let 

l^niV) = P{{iOn,u;n+i,-..)(^V}. (3.28) 



Now by definitions ( p. IS]) , ( |3.17 ), and ( 3.21 ) along with the remark at the end of Section 



3.4, there is some constant < a < 1 such that for every sufficiently large m, 

m 

P{ujj = Xj for each 1 < j < m} = C7([| 6^J(1 + Oia"")) 

for every x G S and some constant C that may be bounded by < ci < C < C2 < cxd with Cj 
independent of x and m. If we define the function : S ^ R to be ip{. . . xiX2 ■ ■ ■) = logftxi, 
then for every x S S, 

P{ujj = Xj for each 1 < j < m} = Cexp{S'„v'(x)}(l + 0(a™)). (3.29) 

Note that this function ip is different from the functions ipp defined in ( p. 17] ) above. Since 
if < 0, there exists a unique 5 > such that the pressure P{6^) = 0. Also, by (^), there 
are some constants < Ci < C2 < oo and a unique ir-invariant probability measure //(^ on 
the Borel sets of S such that 

Ci < ^y(r(a^l^2...Xn^)) ^ 

for every j in the cylinder set r(a;ia;2 . . . Xm)- 

Proposition 9 For A G 7^2 H /C, and every n > 1, the measure fin is absolutely continuous 
with respect to p^ and Pn ^ ^^ip as n ^ 00. Furthermore, P{p) = 0. 



Proof. Our proof is much the same as the one for Theorem 5.4, [11|. For each m > 1, 



the sum over all X1X2 ■ ■ ■ Xm in Gm of the probabilities P{u:j = xj for each 1 < j < m} is 
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equal 1. Since C in ( p.29| ) is positive and finite, by ( |3.30| ), probability measure and 

the pressure P{ip) = 0. Since ( p.29 ) and ( 3.30| ) hold for all cylinder sets and these generate 



the Borel u-algebra, it follows that //i <C fj-^ and that the Radon-Nikodym derivative 
h = {dfii/dfi^) is bounded away from and cxd. Next consider the restrictions of /ii and 
fi^ to the (7-algebra J^n generated by the coordinate functions xj, j > n. Since the tail field 
•^oo = r\n>i -^n is — 1 under fi^ ( fi^ is mixing, see e.g. |^]), we obtain 



Hence, it follows that /i„ <C /i,^ and /i„ ^ /i^. □ 

Thus, the stochastic process u;i,u;2,... is asymptotically stationary, that is, the joint 
distribution of w„,cj„_|_i, . . . converges to that one of a stationary process as n — > cxd. The 
limiting process is a Gibbs state, hence, isomorphic to a Bernoulli shift, but in general non- 
Mar kovian. From ( 3.29| ), we see that each stationary distribution decays exponentially in 



the distance m from the root. The rate at which the measure of a cylinder set T{x\X2 ■ ■ ■ Xm) 
decays depends on the relative frequency of the generators in the reduced word xiX2 ■ ■ ■ Xm 
and the order of the letters. It is easy to show that the distribution is spherically symmetric 
if and only if the contact process is isotropic. 

3.6 Proof that rj < 1 away from K,^ 

The crux of this section consists of finding bounds for the collection of "reinfection" prob- 
abilities. In this section, we restrict ourselves to A G int(/Co)- Let Rxk+iXk+2---^k+m denote 
the event that there is an infection trail that begins at xiX2 ■ ■ ■ Xk at some time, reaches 
vertex xiX2 ■ ■ ■ Xk+m at some later time, and returns to xiX2 ■ ■ - x^ in finite time. Thus, 
Rxix2...xm denotes the event that there is an infection trail that begins at the root at time 0, 
reaches vertex xiX2 ■ ■ ■ Xm at some time, and ever returns to the root from xiX2 ■ ■ ■ Xm- First, 
we note the following about the "infection probabilities" bj. Thanks to constant updating 
while the infection moves off to the boundary of the tree, the hj in the limit approximate 
the actual infection probabilities to any desired degree. Specifically, for large k and fixed 
segment xiX2 ■ ■ ■ Xk-i of some infinite word in S, the probability that xiX2 ■ ■ ■ x^-i • . . x„ 
is ever infected, given xiX2 . • • a^fc-i is infected, is approximated by n^=fc as n ^ oo, 
where each b^^ depends on the initial infinite string . . . xiX2 ■ ■ ■ x^-i- The products nj'=fc ^xj 
approximate the conditional probabilities UxiX2...x„ /uxiX2...Xk_i in the sense that, for large 



k, the sum on the righthand side of (3.15) is the same as the very same sum with 0^=^ bxj 
replacing UxiX2...xn/uxiX2...Xk^i- In addition, for two distinct vertices x = ...xiX2--- and 
y = . . . yi?/2 • • • in 5] and for all sufficiently large m, the probability of the event that 
xiX2 ■ . . Xm or yiy2 . . .ym ever infected is approximately equal YYJL^ ~^ H^Li by^ , where 
the bxj depend on x € S and the by^ depend on y £ T,. 

If we recall the symmetry assumption of the contact process, ( |3.15| ) for p = 2, definition 
( 3.22| ) of the bj, and the facts that the bj are based on the information contained along 



an infinite path, whereas we view a finite piece of the path, a similar reasoning as above 
applies to the probabilities P{Rxk+iXk+2-Xk+m}- Thus, the probabilities P{Rxk+iXk+2...Xk+m} 
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of reinfection are approximated by the products JJ^^k+i^xj (depending on xi...Xk) in 
such a way that the entire collection of probabilities P{Rx\X2---xm 

} is approximated to any 

desired degree in the limit. Hence, from ( |3.18| ), for each A such that r„(A) < 2, for every 
X = . . . X1X2 . . . G S, and all integers m,k > 0, 

P{ U Rx,+,...x,+J = CeKp{SmMx)}{l + Oia'')) (3.31) 

Xk + l---Xk + m,^Om 

for some constants < a < 1 and < C < 00. We point out that an alternative approach 
(not pursued here), which leads to the same estimates, is based on approximating the set 
of vertices ever to be infected by labelled Galton- Watson trees. 

Proposition 10 For each A G int(lZ2 H /C) such that ru{/^) < 2, there exist constants 
< Ci,C < 00 such that for every e > 0, sufficiently small C > 0, and fixed sufficiently 
large k, 

00 

Ci exp{n^(^2 - e)} < -P{ U IJ Rxi...xm} < Cen^; exp{n^(^2 + (3-32) 

m=k Xl...Xm&Gm 

where n^^ = log C • 

Proof. First observe that, by subadditivity, P{Rxi...xm} < P{Rxk+iXk+2-Xk+,r,} (because 
Uxi...x„^ < Uxi...Xk+m,/uxi...Xk)- Fix such A, some e > 0, and some sufficiently small > 0. 
Recall that every generic x £ ^2(0 satisfies 

exp{n^(^2 - e)} < exp{S'„^(/?2(x)} < exp{n(;(^2 + (3-33) 

Next fix some sufficiently large integer k. Then by the previous remarks, by Proposition |^, 
by ( lOll ), and ( p3D , we obtain 



^{U U ™} < ^{U U ^-.+1..-.+^} 

m=k Xl...Xm&Gm rn=l Xk + l-.-Xk + m^Gm 

00 

- '^2 X! exp{Sm^P2{x)} 

m=\ one xx...xm^A\{C,) 

< C2C exp{Sm^2{x)} 

m:|l— r?i/n^|<e One xiX2...Xm(:A2{C) 

< C2C (2en^) exp{n^(^2 + e)} 

for some constants < C2,c < cxd. 

The reverse direction is straightforward. Indeed, by Proposition ^, by (3.31) and ( 3.33| ), 



p{[j U ^-1..-™} > ^{ U ^X™+l..-2^} 

m=k Xl...Xm&Gm XiX2...Xm'^A2{0 

> Ci Y exp{5mV'2(2;)} 

one xiX2...xm&A^{0 

> Ci exp{n^(^2 - e)} 
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for some constant < Ci < oo. This completes the proof of ( |3.32| ). □ 

Time-dependent Infection Probabilities. Next, we turn to discuss the time-dependent 
infection probabihties 

u^^tW = P{x G At} (3.34) 

for X € Q and every real t > 0. Again, the strong Markov property and the monotonicity 
and homogeneity properties of the process imply that Uxy,s+tW > Ux^sWuy^tW for all 
s,t > and x,y £ G such that \xy\ = \x\ + \y\. Clearly, Ux^tW < Ux{X)- A subadditivity 
argument shows that, for every x G Qk and yn = xx . . . x G Qnk-, the limit 

Ymi^Uy^,kns{Xf''' = Ux,s = Ux,si\) (3.35) 

exists, and that Uy^^^nsW < f^:,s(A)'^ for all integers n > 0, in particular, for x = i £ ^, 
we have Uy^^nsW < Ui^sW^- Obviously, since Ux^tW < Ux{X), we have 

Ux,sW < /3x(A) (3.36) 

for every x G Q and each s > 0. Observe that the functions Ux,s{-) share the same basic 
properties with the functions Ux{-), which were essential to prove their geometric decay 
in the distance from the root, for instance, the Holder conditions, and the convergence 
of their sums. Therefore, a parallel analysis may be carried out to define potential func- 
tions and the corresponding Gibbs states as described for the u^(A). To save space, we 
shall omit the details with the exception of a key ingredient to the proof of 77 < 1, which 
is the scaling properties between the space variable x G Q and the time variable s in 
C/x,s(A). To wit, let ^p^six) = ^p^s;x{x) be the analogue of ^Pp{x) = ipp-x{x), where in defi- 
nition ( |3.15| ), on the righthand side, the ratio {uxj^x2...xnW /uxix2.--x,^_i{X))^ is replaced by 

{UxiX2...x„,ns{X) /UxiX2...x^_i{X))^ , that is, 

{Hp4n;x,X2...x,.,;X)),,= ^ [^ ^x,......„,ns(A) y ^33^^ 

UxiX2...Xk-i V-^j 



for every i,j G A. Then the same machinery leads to the potential function ^p s and its 
average with respect to the probability measure ^s^p s > 

^P,s = J <^p,sdns^^ ^. (3.38) 

For every s > 0, let Rxix2-.xmi''^^) denote the event that there is an infection trail that 
begins at the root at time 0, reaches vertex X1X2 ■ ■ ■ Xm at some time, returns to the root 
from X1X2 . . . Xm prior to time ms and stays without recovery mark * up to time ms. By 
the same token as ( |3.32 ) was derived, for every e > 0, sufficiently small C > 0? and fixed 



sufficiently large k, we obtain 

00 

C3 exp{n(;(¥2,s - e)} < P{ y |J Rxr...x^ {ms)} < aenc^ exp{n^(¥2,5 + e)} (3.39) 

m=k Xl...Xm&Qm, 
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for some constants < C3, C4 < cxd and where = n^^^ = log C/<J>2,s- Therefore by definition, 
for every K > 0, 

00 

C3 exp{n^($2,/fs - e)} < -P{ U IJ Rxi...xmi^sK)} < C4en^ ex.p{n^{^2,Ks + e)}- 

m=k Xl...Xm&Qm 

Let denote the time for which the expected number of infection trails returning to the 
root from distance from the root and staying without recovery mark * up to time is 
maximal. Then in view of n^-K = Kn^ and by the additivity of the functional Sn^2,s^ 

CO 

C3exp{n^K($2,tj - e)} < -P{ IJ |J Rx^...x„,imt(;K)} < c^enc^ exp{n^K($2,i^ + e)}- 

m=k X-i_...Xm&Qm 

Combining the last two displays brings K^2,tQ = ^2,KtQ - In other words, distance and time 
are scaled by the same factor in the mean functionals $2,t^) that is, there is some constant 
< A(A) < 00 such that, for each C > 0, 

tc_ = A(A)n(^. (3.40) 

Observe that, due to the form of ( p. 27 ), for each s > 0, 

nQ,s^2,s = riQTp^ = logC- (3.41) 



Proposition 11 For each A in the interior of IC D 7^2 with ru(A) < 2, there is some 
< A = A(A) < 00 such that rj^ = ri{\)^ = exp^2 A) particular, 

exp^2 A < ^ 'i'f o^nd only if 77(A) < 1. 

Furthermore, for each A in the interior of IC with r„(A) < 2, we have r]{\) < 1. 



Proof. First recall that for A in the interior int(7^2 H IC) with r„(A) < 2, we have ip2 < 0. 
Since the function i](-) is nondecreasing in each variable Xj, it remains to be shown that 
77(A)^ = exp^2;A for A G int(7?.2 n fC). 

If we fix some e > 0, some sufficiently small C > 0, and some sufficiently large k, by 
( 3.40| ) and Proposition p^, 

rj = l\m{P{l G At})'/' = l^(p{\J [j Rx,...x^{t)}] 

\ m=k XiX2...Xm<^Gm I 

< lim {{CtuA exp{nA(^2 + e)}]^^*^ 

= ^lim^[(C7enc) exp{nc(^2 + e)}]'/^"« 

= exp{^2 + e}^^^ 
= exp{(^2 + 
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where the employed constant C is positive and finite. 

Next, it is not difficult to see that Proposition |l^ has a time-dependent version (the 
proof being the same) , namely, 



Ci exp{n(^^^($2,s - e)] < P{\j |J RxL^xmij^s)] < Cen^,^ exp{n^,s(^2,s + e)}, 

m=k Xi...Xm&Gm 

for every s > 0. Then the lower inequality follows from similar arguments in combination 
with the product structure of the probabilities Ojifc+i > which the approximation 
( 3.32| ) is based, and again the equality n(^^s^2,s = ^^(^2) that is, 

limiP{l G At})'/' = }}^[P{[j U ^-1.-^(0} ' 

\ m=k Xl...Xm&Qm 

> lim [Ci exp{n^,tA^2,t, - e) 
= hm [C7iexp{nc((?2-e7)}]'/*^ 
= lim [Ciexp{n^(^2-e7)}]i/'^"c 
= exp{(^2 - £7) /A}, 

where 7 = ^2/^2,*^) the constant Ci is positive and finite and k was chosen suitably large. 
Since e > and C > were both arbitrary and Tp2 < 0, the inequality 7] < 1 follows together 
with all other claims. This finishes our proof. □ 



4 First-Passage and Backscatter Matrices 



Recall from ( 3.21| ) that, for fixed x, p, and A, -Bp(A) is a matrix with entries bj{\Y , where the 



bj{X)^ are associated with (p2-\- Now, if we let 6^ = bj{X)'^ be the matrix entries associated 
with the function ^2 A! where ^2-A defined in ( |3.25D for p = 2, then, for any positive 
number p, define the matrix Mp = Mp{\) to be the 2d x 2d matrix, indexed by elements of 
A, whose entries are given by 

{Mp{\))ij = bjiXr ifj/r\ (4.1) 
= ifj = i~^ 

We will refer to Mi (A) and M2(A) as the first-passage matrix and the backscatter matrix, 
respectively. These two matrices will play a distinguished role, for instance, they are im- 
portant in the proofs of Theorems ^ and ^, stated in the Introduction. Note that M2 
coincides with the matrix in ( |1.14| ). Since by construction, for p > and A G /Cq, Mp{X) 
is an aperiodic, irreducible, and nonnegative matrix, thus, a Perron-Frobenius matrix, the 
Perron-Frobenius theorem lets us conclude that Mp(A) has a largest positive eigenvalue 
9{p; A). Henceforth, we will recur to the shorthands 6p for Op{X) = 6{p; A), in particular, we 
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write 



01 = ei{x) = e{i;X) = e{x) (4.2) 

02 = 92{X)=e{2;X). 

Note that 82 = exp^2 that the eigenvalue 9 in Theorem ^ coincides with 9i. In addition, 
it is not difficult to show that 9i = exp^^^ for 1 > r^iX) and 9p = expTp^ for p > r„(A). By- 
definitions (|3.21| ) and ( [4.1] ) and by Proposition 0, the entries of Mp satisfy the equation 

E ^ = 1- (4-3) 

4.1 Strict Monotonicity of the Lead Eigenvalues 

Proposition 12 Let A G /Co and let X he an infection parameter in the direction of decrease 
for A, as defined in Section 2.5. Then for p > ru{X) and every s > 0, 

^p;A < ^p;A, (4.4) 

^p,.;A < ^P.^;A- (4.5) 

In particular, 

9i{X) < 9i{X), (4.6) 

02(A) < 02(A). 

Moreover for each s > 0, each X £ fCo, and t > p > ru{X), we have 

'ft;\ < 'fp-x, (4.7) 

thus, 9t{X) < 9p{X), in particular, 02(A) < 0i(A). 



Proof. Only the proof of ( [4.4| ) is carried out. The proof of ( |4.5| ) runs in parallel whereas 
the second line of (^^) will then follow immediately from the definition of 02 and the first 
line of ( [4.6| ) will follow by recalling the definitions of the bj and 0i. The verification of ( [4.7D 
is deferred to the end of this proof. Fix A G /Co and p > ru{X). Write 



hpix) = hp{x,x,k) = { I I } (4.8) 



n — k — 1 



l/{n-k) 

ri 

UxiX2...x^_l (A) 



/ '^XlX2...Xn{^) ^ 



for all integers n — 1>A;>0,xGS, and A G /Co- Let A be an infection parameter in the 
direction of decrease for A. From ( ^.17 ), we see that there is some < C(2;, A) < 00 so that 



for all sufficiently large n, 

W{xYHp{n; X1X2 ■ ■ ■ Xk-i, A)F(cj""''x) = C{x, X) exp{S'„„fe(/9p;A(2;)}. 
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It suffices to prove that, for fixed x S S, 



lim hp{X,x,k) < lim hp{X,x,k) (4-9) 

(where the hmits exist in view of our assumption p > r^). The meaning of the limit lim„ ,^^00 
is that we take both variables to 00, yet n much faster than k. Once ( |4.9| ) has been shown 
to hold, it will follow that, in particular, it is valid for all x such that Sn^p-\{x) /n — > '^p-^x 
as n — > 00, consequently, '^p.\ < ^p;A- 

Fix X = . . . X1X2 . . . G S. From (|2.6| ), recall that there is some constant < w < 1 such 
that n^i...x„(A) < itxi...a;„(A)w". Choose n sufficiently large so that 



Now combine these to 

UxiX2--x„{^) Ux^x2...x„{X) Ux^x2-..x„{^) Ux-^x2...Xk_l{^) 



U 



'XIX2...X 



k-lW Ux^X2...Xk-lW Ux^X2...Xn{X) UxiX2...Xk-l{^) 

(A) 

m 

UxiX2...Xk-i \'^) 

_ Ux-^x2-..XnW , ,n/2 
— lU 

UxiX2...Xk-i (A) 



Hence, 



and 



hpCX,x,k) < cj^"/(2('^-'=))/ip(A,x,A;), 
lim hpCX,x,k) < lim x, fc) 

n,fc— >oo n./c— »oo 

< tj^/^ lim hp{X,x,k) 

n,fc— +00 

< lim hp{X,x,k), 

which proves ([4. 4 ) since p > ru{X) was arbitrary. Consequently, Proposition implies that 
the function bj{X), j £ A, may be chosen as strictly decreasing function along directions of 
decrease, and similarly, the entries of Mp(A). Therefore, we also verified the two inequalities 



in (4.6) in view of the definitions of 61 and 62. 

Finally, the first statement in (^]^) (and similarly, its second statement) rests on the 
following observation. By arguments along the lines used in the proof of Proposition |^, 
there is a constant < 7 < 1 so that for each A: > and sufficiently large n, we have 
Uxi...x„W/uxi...Xk-iW < 7"- Hence, 

^ 'U'XiX2...x„{^) Y _ ^ UxiX2...x„{^) y-p ^ Uxj^x2...x„{^) ^p 



UxiX2...X^:_l{^) UxiX2...X^._^i{^) UxiX2...Xk_i{^) 

^ „,n^ ^a:iX2...a:„(A) ^p 
UxiX2...Xk-i (A) 



35 



for some constant < 7 < 1 because t > p. Proceeding then along the above route to prove 



(4.4) provides the desired result and ends our proof. □ 

Observe that this proof also shows that the functions in ( |4.4| )-( p^ are nondecreasing 
in each variable Xj because, clearly, for all A and A in /Co with all infection rates Aj = Aj 
the same except for one Xj < Xj, for sufficiently large n, the ratio {ux^...xnW/'^xi...Xk-iW) ' 
(A)/n (A)) < 1. 

5 Continuity Properties of the Infection Probabilities 

This section is concerned with the properties of continuity and discontinuity of the lead 
eigenvalues of the matrices Mp and the functions f3x and ends with the proofs of the results 
about the weak survival phase (Theorems || through |^ and of Theorem |6[ 

5.1 Continuous Potential Functions and Lead Eigenvalues 

If for every e > and for every subset A C R'^, we let A^^ be the e- neighbourhood of A, i.e. 
consisting of all those z S R'^ that have (Euclidean) distance to A less than e, and let JC^ 
denote the complement of /C in R"^, then define /C/ = /C \ (/C'^)^ to be the set of points in /C 



which are at a distance to /C^ at least e. The next result along with Proposition |19| below 
are crucial moments in proving Theorem |l[ 

Proposition 13 For each p > r„(A), s > 0, and all i,j £ A, the functions 

Xj ^ 6i(A) (5.1) 

Aj ^ ^p=¥p;X 
Xj = ^p,s;X 



are continuous for each X G int{IC), with 6j(A) being defined in The meaning of the 

first statement is that each bi{X) can be chosen as a continuous function in each variable 
Xj . In particular, 

Xj ^ iiiX) (5.2) 
is a continuous function for each X G int{IC n 7^2) with ru{X) < 2. 



Proof. Since the verification is the same for both functions and ^*p,s, we only present 
the proof for the former. As explained in the proof of Proposition we can conclude that 



the property holds for bi{X) as well, which will accomplish the first line of (5T). Also, by 
virtue of Proposition O, claim (|5.2D will follow immediately. 



Let p > r„(A). Let e > 0. Observe that for A G fCj, each A^ > 0. For each T > 0, define 

= u^{X) = P{x G At for some t G [0, T]}, 
= uf'°^)(A) = P{x G for some i > T}. 
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Denote hp{X) = hp{X,x, k) as it has already been defined in ( [4.8|) and define 

9pW = 9p{\ X, k) = hp{X, X, k^-^. 

Write K^{\,x,k) and gJ{X,x,k), respectively, for the functions that result when all Ux 
are replaced by in the expressions for hp{X,x,k) and gp{X,x,k), respectively. With 
these agreements along with the observation relying on ( p.lTl) , as described in the proof 
of Proposition 12, in order to prove that the function Xj Tpp = ^p-x is continuous, it is 
sufficient to show that, for fixed x £ T,, ior p > r^iX), and each j G A, the function 

Xj ^ lim hp{X,x,k) = lim lim h^{X,x,k) 

>oo rijfc^ooT— >oo ^ 

is continuous for each A G /C/. Here, lim„^fc^oo means that n and k are taken to oo, n much 
faster than k. Since e > is arbitrary and in view of Proposition 0, our first two statements 
in ( |5.l[) will then follow. 

Key ingredients to the proof are the following three items that we will verify below: 

(A) For each j G A, the function Xj /ij (A,x, k) is continuous for A S Kj and fixed T. 

(B) For fixed A, the convergence of /iJ(A, x, k) to /ip(A, x, k), as T ^ oo, can be controlled 
uniformly for X £ ICj. 

(C) For fixed A, the convergence of hp{X,x,k) to its limit, as n,A; ^ oo, happens in a 
uniform fashion as well for A G /C/. 

These three claims have some consequences. From the first and second claim, it will follow 
that the function Xj — > hp{X,x, k) is continuous. The third claim will imply that the limit 
Xj — > lim„ fc^oo hp{X, X, k) is a continuous function for A G /C/. 

Indeed, the first claim is straightforward since it is obvious that is a continuous 

function in each Xj and that, because each A^ > 0, we have li^(-) > 0. Hence, the function 
hJ{X,x,k) inherits these properties from the 

To show the second claim, note that by monotonicity, for each A G /C/, 

hp{X,x,k) < sup max hp{X,x,k) (5.3) 
gp{X,x,k) < sup max (7p(A, X, fc). 

Next, since u'^-^^x2...Xk-lW — ''^xiX2...x^_i{^), we obtain by Taylor expansion of (1— _fiii±n.y^ 

^\ i\ Tf\ i\ ^ k'^x\X2---Xn^'^^^ ~ i^^i-x-iX^-.-x^X^^^ 

0<gp{X,x,k)-gp{X,x,k) < }^ 

Xfe+ix,+2...x„_iGS„_fc_i {Ux^X2...x,_, lAJ j 

_ {UxiX2...Xn i^))^ 1^ UxiX2^..Xn (^) 

(T,oo) 

/ UxiX2...Xn(^) -.^ 

UxiX2...x„ {'^J 
< CTgp{X,X,k), 
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the reason for the last hne being the convergence of the Taylor series for sufficiently large 
T. Clearly, Ct ^ as T ^ oo. Continuing and using Taylor expansion again brings 



0</ip(A,x,A;)-/iJ(A,x,fc) = (5p(A,x,A;))i/("-'=)-(5j(A,x,A;))i/("-'=) 

gp{X,x,k) 



(5p(A,x,A;))i/("-^){ 



1 gpi\x,k) - gJ{X,x,k) 



n-k gp{\,x,k) 



_^ , 1 gp{\x,k) - gJ{X,x,k) 
n-k gp{X,x,k) 



< hp{X,X,k){^—CT + o{^—CT)}, 

n — k n — k 



where in the last line the upper bound of the preceding calculation provided. In view of 
(|5.3|), the obtained upper bound in the last display is uniform for X G fCj. This completes 
the proof of the second claim. 

Finally, to see the validity of the third claim, recall from (3.18) that there is some 
constant < a < 1 so that 

gp{X,x,k) = C exp{Sn-k^p{x)}{l + 0{a'')) 

for every x = . . .xiX2 ... € S, where the implicit bound in the O(-) term is uniform in 
X and the constant C may be bounded by Ci < C < C2 with Q independent of x and 
n — k. Hence, Sn-k^p-\{x) / {n — A;) — > ipp.\{x) as n,k ^ 00. Considering those x such that 
Sn^p;\{x) /n —>■ ^p-\ as n,k 00, and combining the last display with ( |5.3D completes our 
proof. □ 



Observe that for the proof of the statements in (5.1), it was not necessary to use rj < 1. 
The proof shows that, for p > r„(A), the function Xj bi{X) is continuous as well for each 
A G int(/C) and i,j G A. 

Corollary 1 For each j G A, the function Xj — > VuiX), defined in ( 3.1(\ ), is continuous for 
each X £ int{IC). 



Proof. First recall that to each p > there is associated a matrix Mp{X) whose entries 
bf.{X) may be chosen as continuous functions in each variable Xj by Proposition 13. In that 
case, all matrix entries bk{X) are less than one. From these facts, it easily follows that p 
can be varied continuously so that the norm, say, of Mp{-) stays constant as Xj is changed 
to A'- in a neighbourhood of Xj, and thus, ru{X) must be continuous in Xj. □ 
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5.2 Mean Offspring Number /i^ of the Galton- Watson Trees 

Recall Wx = PiV^} from ( |2.1| ). Since Wx{X) < Ux{X), it follows that the collection of 
functions Wx satisfy a Holder condition as do the functions Ux- Consequently, with each 
function Wx{X) there is associated a well-defined potential function constructed by 

the same means as described to obtain {pp-x{x) corresponding to Ux- Next we show that the 
two potential functions (x) and ^p{x) coincide for every p > ru{X)- 

Proposition 14 For every A € /Co and p > r„(A), the two potential functions and ipp 
coincide, that is, for every x G S, 

(^-(x) = (5.4) 

in particular, 

^P=^p- (5.5) 



Proof. Let A € /Co- Ii^ order to facilitate the presentation, we may carry out the proof for 

fc "^X-i_...Xn 



p = l> r„(A) and look to bound Exfc+i...x„eg„_fc Wx^...x^ instead of Y^Xk+i.-xr^ag. ■ 



Minor modifications establish claim (|5.4| ) for arbitrary p > ru(A). By definitions ( 3.15| ) and 
( 3.17] ), we collect for all integers k < n and each x S S, 



Xk + i...X„&gn-k Xk^i...Xn&Qn-k 

< {max Uy) 



< (maxujy) Cexp{S'n_fc(yfi(x)} 



for some finite constant C. Therefore, for fixed k and each x G E, we obtain 

limsup( Y Wx^...x„y^'' < limsup„^o^ exp{5„_fc99i(x)/n}. (5.6) 



n—^co ^„ 

Xi;^l...Xn&yn-k 



It remains to be shown that, for all sufficiently large k and x G S, 

liminf( ?^xi...x„)^^"' > liminf exp{5„_A;V?i(x)/n}. 

T) — ^nn ' . r) — ^rvi 



n— >oo 

a;fc+i...x,iGe„_fc 



This will imply claim ( |5.4| ) because lim„_»oo 'S'„_fc(/?i(x)/n = (/Ji(x) exists for 1 > ru(A). 

We shall mimic the construction of the proof of Proposition 1 in ||T^. For each x G Q.^, 
let V[x,k} denote the probability that there is an infection trail from the root 1 to x that 
remains within distance k of the geodesic segment from 1 to x. Clearly, v^x,k} — and, as 
k ^ oo, V{x,k} T Ux, that is, limfc^oo V{x,k} = Ux- 

Thus, by definitions ( |3.15p , ( p.l7| ), and ( p.21D , if we write na;,.,,^;^ = Ux^...Xkiuxi...x„,/uxi...Xk), 



then for any e > 0, for all m and k sufficiently large. 



V{x,k} > - ^)uxi...xk n (5.7) 

j=A:+l 
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where each bxj = bx^ (A, xi . . . Xk) depends on xi . . . Xk and m. 

For every vertex z G Qk at distance k from the root 1, define to be the probabihty 
that there is a direct infection trail from the root 1 to z, that is, a trail which follows the 
geodesic segment from 1 to z. It is apparent that > 0. 

Now choose x = xiX2 . . . x„ G ^„ so that n = 2k + j + Nm for some < j < m — 1 and 
let yo,yi, ■ ■ ■ ,yN be the vertices on the geodesic segment from the root 1 to x such that 
Vi G Gk+im- Suppose that all of the following events take place: (A) There is a direct infection 
trail from 1 to yo that reaches yo at a stopping time Sq, (B) for each i = 0,1, . . . , N — 1, 
there is an infection trail from yi to yi+i that begins at time Si and ends at time 5'j+i, 
which remains within distance k of the geodesic segment from y^ to yi+i, (C) there is a 
direct infection trail from y^ to x, beginning at time Sn- (If we let Si be the first time 
after Si-i that such a path arrives at yi, then the random times Si are stopping times.) 
Observe that concatenating the infection trails (A), (B), and (C) constitutes a downward 
infection trail from the root 1 to the vertex x. Because the Si are stopping times and the 
events (A), (B), and (C) occur on nonoverlapping parts of the percolation structure, by the 
strong Markov property and the monotonicity properties, we have 



N-l 



i=0 

Notice that yo{yQ^yi){yi^y2) ■ ■ ■ (yivl^iJ/Ar) = yN- Hence by (|5.7|) , we have shown that 

n—k—j 



Wx>ay,,a -i{l-e)^{ux,...x^)^ H ^■ 



3=k+l 



If we write ai = min^^^g^ Uy, 02 = min^gg^,^^, ay, and = min^.gg^ Ux, by ( 3.17 ), we obtain 



n—k—j 

€Q„-k Xfc+i...3::„e5„„fc jf=fc+l 

> 0102(1 - e)^(u*)^ C exp{Sn-2k-j^iix)} 

for some positive constant C. Note that the employed constant C (as well as C earlier) 
takes care of the discrepancy in the index sets of the summations. Then 

liminf ( V 'Wxi...x„)^^" > [(1 - e)n*]^/"^ liminf exp{Sn-2k-jV>i{x)/n}. 

As e > is arbitrary and as m can be chosen sufficiently large for suitable k, the desired 
result (5.4) now follows. Finally, (^.5|) is an immediate consequence of definition (^.25). □ 



Lemma 5 As n — > 00, 

( n.(A))i/"^exp(^i;, = 0i(A). (5.8) 
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Proof. To simplify our exposition, we will present the proof in terms of the Ux and V3i;A, 
rather than lifting the calculation to and f2:\- However, observe that the two calculations 
are related via the functions bj. For every integer A; > and a; G S, by subadditivity and 

"^^^^^ - ~^ or 

< C exp{Sn^i;x{x)} 
for some finite constant C. Hence, for every x £ T,, 

limsup( y2 ^x'('^))^^" < limsup exp{S'„i^i;A(a^)/''^}- 

X&Qn 

Since this inequality holds for every x G S, it must hold for those x £ T, with Sn(pi{x)/n 
Tpi as n ^ oo. Thus, letting /c — > oo, it is easily derived that 

limsup( V n^(A))^/" < exp^^.;^. 

n-->oo 

xey„ 

To show the reverse direction, fix e > 0, and define 

= {xiX2 . . .Xfc e ^fc : UxiX2...xA^) >ex.p{{lp^.x- e){n- k)}}. 

Xk+l---Xn€Gn-k 

It is obvious that for sufficiently large n and fixed sufficiently large k, the set Qs is nonempty. 
Therefore, 

J2 '"^(-^) ^ J2 Ux^X2...x„W 

X&Gn X\...Xk&Qe 

> exp{(^i;,-e)(n-A;)}. 



l^igf ( E ^xWY'"" > - e}. 

xeGn 



Whence, 



As e > was arbitrary, this verifies claim ( |5.8| ). □ 

Corollary 2 // we recall Cn from Section 2. 2, then 

hm ( Wx{X)f''' = lim ( 5] w;.(A))i/" = exp^^.;, = OiiX). (5.9) 

X&Cn XI^Gn 

Proof. Again, our proof is in terms of the Ux and ipi-\ (see the remark at the outset of the 
proof of Lemma ^). In view of Proposition |l^ and Lemma |5|, the leftmost equality in ( ^.9[ ) 
remains to be verified. This may be accomplished by showing that 



hm ( Y ^xf" = lim ( uxf'^. 

XGCn x£Qn 
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But in light of 



J2u, = l'Hi{n;l;-)l, (5.10) 
^ = u*Fi(n;l;-)l, 

where 1 is the vector ah of whose entries are 1 and Ua is the vector with entry in the 
slot and all other entries 1. Since Hi{n;l;-) is a Perron-Frobenius matrix with lead 
eigenvalue exp Snfi{^), taking the n-th root and n ^ co on both sides in both lines of 
( |5.10| ) in combination with Lemma |5| proves claim ( |5.9| ) . □ 

Corollary 3 The mean offspring numbers fik for the Galton- Watson trees Tk, as defined in 
Section 2.4, satisfy 

lim nl^'' = exp^i = 9i{X). (5.11) 



Proof. The mean offspring number /x^ for the Galton- Watson tree is, by construction, 



This sum differs from the one on the lefthand side of the left equality in (^]^) in that the 
smaller index set £^ replaces C^, thus, the vertices x have word representations ending in 
the letter a. Consider a vertex x £ Ck- Clearly, there is positive probability p, independent 
of X, that a vertex y £ ^ '^(^) infected within two time units. This brings 



Corollary |2| thus finishes the proof of ( ^.11 ). □ 



For all integers r,n > 1, define Zn{r) = Zn = |Ki('")| to be the cardinality of the n-th 
generation of r^. For n = 1 for instance, EZi{r) = X^xg^* 

Lemma 6 Suppose that Oi{\) > 1. Then for any 1 < 7* < 6'i(A), we have 

liminf P{Zi(2'=) > {-i,f] = /> > 0. 



Proof. This proof is essentially the same as the one in [^] (Corollary 3). We shall put it 
in our context. Let x £ ^2*=- Fix some 1 < m < k, and let = 1, xi, . . . , x be the vertices 
along the geodesic segment from 1 to 2; at distance i2™ for i = 0, 1, . . . , 2^~"^. If for every i, 
there is a downward infection trail from Xi to rcj+i that begins at the time of termination 
of then there is a downward infection trail from 1 to x, thus, x G ^1(2^^). Therefore, 
X e V2k-m{2"') implies x £ Vi{2''). Hence, Zi{2'') > ^2^-^(2^"). By construction, ^(2"") is 
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a Galton- Watson process with mean offspring number EZn{2"^) = X)ze£* m ^ standard 
result from the theory of Galton- Watson processes tells us that 

= z 

exists and, because the offspring distribution has finite support, Z > almost surely on 
the event of nonextinction (see Theorem 2, Section 6). This means that P{Z„(2™') > 
((7*)^'")" eventually} > 0, which together with Zi{2^) > Z2fc-m(2™) concludes the proof. 

□ 



Proof that < 1 at the Transition to Survival. Consider the set of all vertices 
in Qm ever to be infected. In this paragraph and again when we address the behaviour of the 
contact process at the first phase transition, we will exploit the fact that, on homogeneous 
trees of degree larger than 2, whenever £"13^^! — > oo as m — > oo, a (labelled) Galton- Watson 
tree may be embedded in the set of vertices ever to be infected, which gives rise to a 
supercritical Galton- Watson process. This Galton- Watson process grows without bound 
with positive probability, which implies that, with positive probability, \Af \ — > oo. 

Clearly, 

E\ym\ = E ^ jg infected} 

~ $Z ^'^{x is ever infected} 
= 

where /|,} denotes the indicator function. In view of Lemma ||, we know that for A TZi, 
the lead eigenvalue 6i{X) > 1, because ^i(A) < 1 would imply that J2m=iE\ym\ < K 
for some constant K < co, and thus, the contact process would not survive with positive 
probability, which contradicts our assumption that A T^i. As a consequence, r„(A) > 1 
for A E int(7^2U7^3). 



Lemma 7 For A G 7^l n 7^^ n int{lC), 



Oi{\) < 1, (5.12) 
r«(A) < L 



Proof. Throughout the proof, let A G int(/C). It is enough to prove the first claim. For 
d = 1, claim (|5.12| ) is obvious because 1 > /3(A) = ^i(A). Thus, let d > 1. We suppose that 
^i(A) > 1 for A G 7^1 n T^f and proceed by contraposition. By Proposition 13, for each 
j £ A, the function \j — > ^i(A) is continuous. Furthermore, by Proposition 12, for any A* 
in a direction of decrease for A, we collect 6'i(A*) < Oi{X). By a combination of these two 
properties, because 9i{X) > 1, there exists a A* in the interior of TZi such that 0i(A*) > 1. 
Thus, for this A^,, -E|3^m| — > oo as m — > oo. Now, by Lemma |6|, for sufficiently large k (as k 
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runs through powers of 2), the Galton- Watson process {Zn{k)}n>o is supercritical and the 
corresponding Galton- Watson tree has positive probability to grow to infinity. Because the 
Galton- Watson process is dominated by the number of vertices ever to be infected, it follows 
that, with positive probability, the contact process survives. However, this contradicts our 
assumption that A* G int(7^i) C TZi. Hence, it must be the case that 9i{X) < 1. □ 



5.3 Asymptotics of the and Related Functions 

Inspired by some ideas in [^], we collect some bounds between various key functions that 
are equationally interrelated. In fact, our Lemmata |8|, ^ Proposition |l^ and their proofs 



are adjusted versions of the reasoning in |26|. We will prove the results in this section. 



Section 5.4 and in Propositions and ^ under the following Standing Hypothesis, which 
will shortly turn out to be extra (for more on this, see Corollary ^). 

Hypothesis I. Assume that, for every e > 0, the exponent r„(A) < 2 for each A G 
/Cn (/C^)e. 

A consequence of Hypothesis I is that, for every e > 0, we have r„(A) < 2 for each 
A G /C/, due to the monotonicity properties of the contact process, where /C/ = /C \ (/C'^)e 
denotes the set of points in IC which are at a distance to IC^ at least e. 

Next define 

/•oo 

Ux = UxW = P{xe At}dt. (5.13) 
Jo 

Obviously, Ux < Ux for every vertex x £ G. 

Lemma 8 Let e > 0. Under Hypothesis I, there exists some constant < c < oo such that 
for every A G /C/ and x G Gn, 

Ux < c{n + l)ux. 



Proof. Let denote the probability for the isotropic contact process with infection pa- 
rameter e > and initial infection at the root vertex that between time and 1 there 
is no recovery mark * at the root 1 and that, for some i £ A, the root vertex infects its 
neighbour i. Thus, = e~^(l — e"*^) > 0. Fix A G /C/ and fix e* > such that each Afc > 
(e* = e/V^d should suffice). Then for every x = xiX2 ■ ■ ■ Xn £ Gn and y = yiy2 ■ ■ - Vt ^ Gt 
such that \xy\ = n + t, it follows that 

P{leAt+n} > (a,J"P{xGAa. (5.14) 

Fix A G /C n {IC^)e in a direction of increase for A (so that each I3j{\) > (3j{X) > 0). Thus, 
^uiX) < 2. Recall that ri(-) is nondecreasing in each variable Xj and from Proposition 11 



that r/(A) < 1 for A G int(/C). Thus, P{1 G At} < r/*(A) < r]\X) < 1 for every t > 0. This 
observation combined with ( |5.14| ) gives 

P{x G At} < a-"r/(A)*+" < a-"ry(A)*. (5.15) 
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Recall that Ux < Ux- Therefore, if K > 0, by ( |5.15| ), 

roo 

UxiX) = / P{x£At}dt 
Jo 



poo 

< Ux{\)Kn + / P{x G At}dt 

JKn 

< Ux[X)Kn + a-"r/(A)^7log(l/77(A)). 

For K sufficiently large, in particular, such that rj{X)^ /a^^ < 1, we obtain 

Ux{X) < 2KnUx{\) < c{n + l)ux{X) 

for some < c < oo and for each n > 0. The case n = is an easy instance of r/(A) < 
r/(A) < 1. This ends our proof. □ 

For any a; G ^, let Xx denote the total number of infection trails (infection arrows in 
the percolation structure without recovery marks) leading to x and let 

h = h{X) = EXi + l 

Ix = Ix{X) = EXx- 
It is apparent that Ux < Ix for every x £ G. 
Lemma 9 For every x,y €z G such that \xy\ = \x\ + \y\, 

Ixy < Ixly (5.16) 



Proof. Clearly, ( 5.16| ) holds for x = 1 or y = 1 (recall that 1 denotes the root). As was done 



in |£6[, we argue by comparing the underlying process to a multitype contact process that 
evolves as follows. Think of an infected vertex as a site hosting a particle. Infected vertices 
can carry infections of types 0, 1, 2, ... , thus, particles can be of type 0, 1, 2, ... , with the 
initial particle at the root being of type 0. Particles of type 1,2,... evolve as independent 
contact processes, independent of particles of type 0, in particular, particles of different 
types can coexist at the same vertex at the same time. Fix a vertex x £ Q. Particles of type 
evolve as a contact process, except for the fact that they cannot infect vertex x. The first 
time a 0-particle attempts to infect x, a particle of type 1 is placed at x. The second time 
a 0-particle attempts to infect vertex x, a particle of type 2 is placed at x. And, so forth. 

This multitype contact process dominates the contact process in the sense that, when 
there is a particle of the contact process at some vertex, then there is at least one particle 
of the multitype contact process at the same vertex. Let X'^ denote the total number of 
attempts by a 0-particle to infect x (thus, the number of types of particles distinct from 
0), and let X'^y denote the total number of attempts by a particle of any type to infect the 
vertex xy, where y is such that \xy\ = \x\ + \y\. Thus, the chain of inequalities 

Ixy < E{X"y) 

= E{E{X"y\X'^)) 
= EiX'^EiXy)) 

= E{X',)E{Xy) 

< E{Xx)E{Xy) = Ixly 
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finishes the proof of ( 5.16| ). □ 

Similarly as for Ux, a subadditivity argument shows that, for any integer k, x £ Qk, and 
each Un = XX . . . X & Qnk, the limit 

hm /V"=/3, = /3,(A) (5.17) 

n — >oo 

exists and that /y„ > /?" for each integer n > 0. 

Recall that an infected vertex x £ Q attempts to infect its nearest neighbour xa at 
the rate Xa for a G A. Therefore, enlarging on equation (3.2) in yields the following 
recurrence relation, 

4 = ^aUxa- (5.18) 



Proposition 15 Let e > and assume Hypothesis I. Then there are some positive finite 
constants ci(e) and C2(e) such that for each X £ fCj, each integer k > 0, and every x £ Qk, 

h{X) < ci{e)ux{X) < C2(e)(n + l)ux{X) < C2(e)(n + 1)4(A). (5.19) 

For each X G int{IC), each x G Qk, o-nd each yn = xx . . . x £ Gnk, we have 

l™,n,„(A)i/- = lim /,„(A)V« = p,{X) = lim u,„(A)i/" = /3,(A). (5.20) 

I L ^ CXJ 7 1 ^ CXj TL * OO 

Furthermore, there are some constants < Cj(e) < oo, z = 1, 2, 3, such that for X € /C/, for 
all integers n,k > 0, every x £ Qk, and yn = xx . . . x G Qnk, 

%„(A)</3,(A)- (5.21) 

%„(A) <C3(e)(n + l)/3,(A)" 
V(A)<C7i(e)(n + l)/3..(A)". 



/?x(A)" ^ 

Ci(e)(n + 1) - 

C2{e) - 

/3x(A)" < 



Proof. Let e > and A G /C/. Again fix some > such that each Aj > (e* = e/\/2d 
will do). Then by (|5.14|) , for every x = xiX2 ■ ■ - Xk £ Qk and any Xk+i G A, we have 



XlX2-..Xk 



/oo 
P{xiX2...Xk G At}dt 

/oo 
a^^P{xiX2 . . . Xk+i G At-i}dt 

roo 

= / a^^P{xiX2 . . .Xk+1 £ At}dt 
Jo 

= O'e* 1*xiX2...a;fc+i (A) 



(Note that xiX2 ■ ■ -Xk+i may be in Qk-i)- This inequality combined with Lemma |^ and 
relation ( 5.1^ ) implies 

4-(A) < ci(e)u^(A) < C2(e)(n + l)u^(A) < C2(e)(n + 1)4.(A) 
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for each x £ Gk and for some constants ci(e),C2(e), both depending on e, as advertized in 
( 5.19| ). Next apply ( 5.19| ) with x = y„, take the n-th root in each expression and take hmits. 
Since e > is arbitrary, we conclude that, for A G int(/C), 

/3x(A) = I3,{X) 

and obtain ( |5.2[1| ). Combining this with ( |1.6| ) and ( 5.17| ) yields all of ( 5.21| ), as desired. □ 



Corollary 4 Let e > 0. Under Hypothesis I, there is some positive finite constant C = C(e) 
such that for each A S /C/, all integers n — 1 > k > 0, and every x = . . . X1X2 . . . E S, 



where 



g' = g'{X,x,k) = Y: ^""^^^ 



g^ = g^{X,x,k) = Y: (5-22) 



Proof. This is an easy exercise thanks to ( [5.19| ). For instance, there is some positive finite 
constant such that for all integers n > k > and every x £ Gn, 

Ck ttxi...Xfc_i(A) Ixi...Xk._i{X) Uxj^...Xk-iW 
Analogously, the ratio for can be squeezed in. □ 

An interesting consequence of Corollary ^ is that under Hypothesis I the results, obtained 
for the functions Ux in Section 3, can be derived for the functions Ix and Ux, more precisely. 
Holder continuous matrices may be constructed associated with Ix and Ux, respectively, as 
well as potential functions ^Pp{x) and f^ix). But in fact, it is an instance of Corollary ^ 
that these two functions are the same and coincide with the function '^p{x) corresponding 
to Ux, that is, for each x G S, 

9^p;a(^)=<a(^)=<^p;a(x). (5.23) 

Thus, the lead eigenvalues of the corresponding matrices Bp coincide for the three functions 
Ux, Ix, and Ux and their matrix entries may be chosen the same. 



5.4 Continuous Limits via Subadditivity 



In the isotropic case, continuity of the function in A on (0, A2) has been proven in |26|. 
An analogue holds in the anisotropic case with an analogous proof. 
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Proposition 16 Under Hypothesis I, for all j G A and x G G, the function Xj PxW is 
continuous for every A G int{K.) and left- continuous for A E R^|_. 

Proof. The proof for the isotropic contact process, given in [^], carries over with minor 
adaptations. Since the reasoning is short, we present the details here. 

It suffices to show that, for all j £ A and x G G, (a) Xj (3x{X) is left-continuous on 
R'^, and (b) Xj —>■ /3x(A) is right-continuous for every A G int(/C). 

To show the first of these claims, note that, if for each x £ Gk, Un = xx . . . x £ Gkn, and 
every T > 0, we set = Uy^{X) = P{yn G At for some t S [0, T]}, then 

Px{X) = supK„(A))1/" = supsup(n^„(A))i/". 

n>l ra>lT>0 

It is obvious that Uy^{-) is a continuous function in each Xj, and thus, Px{-) is lower- 
semi-continuous in each Xj. Since the function l3x{-) is nondecreasing in each Xj, it is left- 
continuous in each Xj. 

To verify the second claim, for each x £ Gk and for i/n = xx . . . x G Gkn, we may define 



u^JX) = J^ P{yneAt}dt. 



Fix e > and e* > such that each Xj > e*. For fixed A, as T — > oo, the convergence of 
Uy^{X) to Uy^{X) can be controlled uniformly for A G /C/, where ICj was described earlier. 
To see this, fix A G /C/ and A € /C n {IC'^)e in a direction of increase for A, and recall 
a, = e-i(l - e-') > and ( p^ . Then 

0<%„(A)-n^JA) = P{y^eAt}dt 

roc 

< (aej-'^" ri(xf+^'^dt. 



Since, by Proposition 11 and Hypothesis I, we have ri{X) < 1 for A G int(/C), the expression 



on the lefthand side of the last display vanishes as T oo. Thus, Uy^{-) converges to Uy„{-) 
uniformly on /C/. It is apparent that, for each j € A, the function Xj —i- Uy^ (A) is continuous, 
and thus, the function Xj Uy^{X) is continuous as well on ICj. From ( |5.19| ) and the left 
inequality in the second line of ( ^.2lD , we collect iC2{e)uy^{X))^^"' > (3x{X) and, also, from 
Proposition |l^, lim^^oo 

(C72(e)nj^„(A))V" = /3^(A). This implies that 
PxiX) = inf{C2{e)uy^{X))'/\ 

n>l 

Hence, the function Xj Px{X) is upper-semi-continuous on ICj. Since the function Xj — > 
f3x{X) is nondecreasing, it is also right-continuous on this region. Since e > is arbitrary, 
the advertized claim follows. □ 
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5.5 Discontinuity 

From recursion ( [5.18D , we collect the following fundamental relation. 

Proposition 17 Let e > 0. If Hypothesis I holds, then for every A G /C/, every i G A, and 
X = x\X2 ■ ■ ■ Xn-ii G Gn, there is some Dx{x) (independent of i) so that 



bi{X) 



\Dx{x)+b,{X)Dx{x) 

jeA-{i} 



(5.24) 



where the bj = bj{X) are defined in ( 3.21 ) and Ci < D\{x) < C2(n + 1) for some positive 
finite constants Ci and C2, independent of x,n, and X. 



Proof. Fix A G /C/. Suppose that Aj > 0, thus, Pi{X) > 0, as otherwise ( ^.24 ) is swiftly 



verified. Assume that x = xiX2---Xn-ii G Gn for some i G A. Identity (|5.18D may be 
rewritten as 

Ix/u^i'-^ = Xi+ XjU^j/u^i-i, (5.25) 

jeA-{i} 

which in light of (|5.1S| ), the definition of the matrix Bi{X), identity ( ^.23 ) and the remark 
following ( 5.23| ) may be restated as 



bi{X)/Dx{x) = Xi + bi{X) J2 ^i^iW' 

jeA-{i} 



where l/ci(e) < Dx{x) < (n + l)c2(e)/ci(e) because, by (|5.19D , ci(e)/((n + l)c2(e)) < 
Ix/ux < ci(e). Since the constants ci(e) and C2(e) can be bounded uniformly on ICj, this 
ends the proof. □ 

Observe that the system of identities ( ^.24 ) may differ for different A since the coefficients 
Dx{x) are not necessarily the same for different A. For fixed A and x G ^, those are 2d 
equations of which the d equations indexed by j G A-^ are redundant since A, = Aj-i and 
bi = Hence, system ( 5.24 ) may be regarded as a system of d equations in 2d variables 
or in d variables. 

Proposition 18 Let e > and assume that Hypothesis I holds. Then for every A G /C/, we 
have l3i{X) = f3j{X) for all i,j €z A if and only if Xi = Xj for all i,j £ A (that is, the contact 
process is isotropic). Equivalently, bi{X) = bj{X) for all i,j £ A if and only if Xi = Xj for 
all i,j G A. 



Proof. It is clear that the statement in terms of the is equivalent to the statement in 
terms of the bi. Thus in this proof, we will restrict our attention to the collection of bi. 

Pick A G /C/ and fix x G ^. Write D = Dx{x) > and Cx = J2jeA'^j- Fii'st, suppose 
that bi{X) 



b(X) for each i G A. Then in view of Proposition [l^, relation ( 5.24 ) reads 



6(A) = XiD - XiDb{Xf + Db{Xf ^ A^ 

jeA 



X,D{l-b{Xf) + Db{XfCx 
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for every i G A. But, because this holds for all i €z A, it must be the case that, for any two 
distinct indices i,j G A, 

XiD{l-b{Xf)=XjD{l-b{Xf). 

However, if for some pair i,j G A, it was true that Aj ^ Xj, then this would mean that 
6(A)^ = 1, thus, 6(A) = 1. But this contradicts our assumption that A G /C. Therefore, we 
conclude that all infection rates Aj are identical, thus, the contact process is isotropic. 

To show the reverse direction, suppose that the contact process is isotropic, that is, 
Aj = A* > for all i £ A. Then equation ( ^.24 ) assembles as 

bi{X) = X^D - X^Dbi{Xf + X^Dbi{X) bj{X). 

jeA 

Equivalently, if we write Cb = J2jeA^jW^ 

fe.(A)^ + b.(A) ^~f^^-^ -l = 0. 

Thus, we find an explicit expression for 6i(A) by means of the quadratic formula, which is 
the same for all i G ^ (because the coefficients of the polynomial are independent of i). 
Hence, 6j(A) = 6(A) for all i G ^, as required. □ 

The following Proposition is crucial to prove Theorems [l|, ^ and |^. 

Proposition 19 There exists a set T>c C dIC such that, for each X G "Dc, the {6j(A) = 
bi{X,x)}i^A satisfy 



Proof. In it was shown that, for nonnegative constants pi < 1, not all zero, with all 
Pi = Pi-i, the system 

Fi{z) = zpi - zpiFf{z) + Fi{z) J2 zpjFj{z) (5.27) 

of 2d equations is satisfied by analytic functions {-Fj(2;)}jg^ within the largest circle not 
containing any singularities, and, that the system has a singularity for z = R, where R is the 
radius of convergence of the power series representation of Fi for each i G A. Furthermore, 
each of the functions Fi{z) is an algebraic function, i.e. it satisfies a polynomial equation 
with coefficients in the ring C[z]. (See Q, Section 3.2, for a discussion on an elimination 
algorithm and resultants to obtain Fj.) But system ( ^.27 ) is the same as system ( 5.24| ), 



where now the 6j(A) play the role of the Fi{z) and each XjDx the one of pjZ. Therefore, 
because there does not exist any continuous choice for the functions 6i(A) at the singularity 
of the Fi{z), it is clear that if the system is singular for the Fi, it must be singular for the 
6j. Because the argument is an important piece in our approach and is short, we will derive 
the equation that characterizes the singularity. 
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Let Vc denote the set of singularities, that is, the set of A such that system ( ^.24 ) is 
singular. If we use the short notation D = Dx, we recah the equation 



hi{\) = \iD + hi{\fD\, + 2bi{\)D ^jbj{X). (5.28) 

jeA+-{i} 

The following derivation is in terms of the hi instead of the Fi to avoid switching notation, 
i.e. set bi{X) = Fi{z). The derivatives are symbolically justified by ( ^.271) . Next we rely on 
the Implicit Function theorem to locate the discontinuity. If {Hi{{hk}k£A+)}i<^A+ denotes 
the difference between the lefthand and righthand sides of ( 5.28| ), then the Jacobian matrix 
of the d X d system of equations is given by {dHi/dbj)iji^_/(^, and whence, by the (complex) 
Implicit Function theorem (see e.g. [|lO|), must be singular at every A G Pc- Since 



1 = dbi/dbi = 2XiDbi + 2 ^ XkDbk (5.29) 

keA+-{i} 

= dbi/dbj = 2XjDbi, U i) 

the Jacobian matrix may be written as I — J(A), where J has entries 

J{X)ij = 2 J2 ^kDbk if j = i, (5.30) 

= 2XjDbi iij^i- 

Since these entries are nonnegative, the spectrum of J(A) is contained in the closed disk 
with radius 7(A), where 7(A) denotes the lead eigenvalue of J(A). In [^, we showed that 
the function 7(-) is analytic, nondecreasing and continuous in its arguments. For A with all 
Aj small, the entries of J(A) are small since they are linear combinations of the Aj, which 
implies that, for A sufficiently close to the origin, 7(A) < 1. Therefore, the A of smallest 
distance to the origin so that I — J(A) is singular must be all those A, where 7(A) = 1. In 
other words, Vc = {X : 7(A) = 1}. 

Next we analyze the corresponding eigenvalue equation J{X)v = v for A G Pc- As a 
consequence of the Perron-Frobenius theorem, the vector v has all entries nonnegative and 
at least one strictly positive. Thus, set s = 2J2jeA+ ^j^'^^j > 0. The eigenvalue equation 
J{X)v = V together with equations (|5.30 ) may be rewritten as 

Vi = Vi[2 J2 XjDbj] + 2bi J2 ^i^^J 

ie^+ jeA+-{i} 

= Vi[2 XjDbj -2XiDbi]+2bi ^ XjDvj 

jeA+ jeA+ 

= 2vi ^ -^j^bj +biS. 

jeA+-{i} 

Multiplying each side by bi and substituting the arrangement of relation (5.24), i.e. 

bi - XiD - XiDbf = 2bi ^ A^- D bj, 

jeA+-{i} 
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gives 



Vibi = Vi [bi - XiD - XiDbf] + bfs 
^ v,DX^ = sb^J{l + b^i). 

Recalling that s = 2^^,^^^^ XjDvj together with the symmetry of the contact process 
immediately provides relation ( ^.26 ). 



It remains to be shown that C dfC. Observe that, in view of the continuity (Propo- 
sition ^) and monotonicity properties of the functions bi, the expression in ( p.26 ) uniquely 



determines the singularity in the reals and the set T>c cannot lie inside the region /C. On 



the one hand, if d = 1, equation ( ^.261) says that bi{X) = Pi{X) = /3(A) = 1, thus, obviously 
A G die. On the other hand, if d > 1, first note that, for A G IC^, there exists some a £ A+ 
with Pa{X) = 1, thus by Lemma ^, with ba{X) = 1. Now choose some A G int(/C). Since 
d > 1, we can assume that 6fc(A), ba{X) > for some a,k £ A. Pick such a and k in A. Sup- 
pose that 26fc(A)^/(l -|- 5fc(A)^) = S > 0. Therefore, due to the monotonicity of the function 
bk{X), for each A G Pc, 



^^^. + 6,(A)2 - 1 + 5,(A)2 



A- 6, 



By symmetry, 6fc(A) < 1, too. As a consequence, since a and k m A were arbitrary, A G /C. 
Combining these gives Vc <Z K. and Vc C dK.. Also note that each ba{-) experiences a jump 
discontinuity for each X £ T>c- This completes our proof. □ 



We point out that, perhaps surprisingly, identity ( p. 26 ) does not involve the Dx that 



appear in relation ( p. 24 ). In other words, the form of equation ( |5.26| ) is independent of the 



D\. Recall our strategy of proof from Section 2.5. We now carry out the last few steps. The 
following observation continues the discussion in the last part of the proof of Proposition 

VcCiC if and only if dlC C JC, 

equivalently, 

Vcd tC if and only if each bi{X) < 1 for X £ Dc- 

We conclude that Vc C IC for d > 1, and, that Vc C /C^ for d = 1 because X £ Vc implies 
A /C. We end this section by collecting some consequences of the preceding result. 



Corollary 5 For X £ Vc 



and for X £ int{IC), 



expv'2;A = (A) = 1, 
VW = 1, 



expy52;A = ^2(A) < 1, 

V{X) < 1. 
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For each j £ A, the functions Xj ^2-\ ^'^^ ~^ vW '^'^^ continuous for A G int{IC), are 
left- continuous for A G /C and d > 1, and, are right- continuous for A G /C^ and d = 1. The 
functions Xj ^'^^ ~^ vW ^'"fi ™^ right- continuous for X £ Dc <ind d > 1 but are 

left- continuous for X £ "Dc CLnd d = \. In addition, all above statements on (dis)- continuity 
are valid for the functions f3j{X) and bj{X) for j £ A, and for ^p-\ for each p > ru{X). 
Furthermore, 

ruiX) < 2 for X£ int{IC), 

ru{X) < 2 forX£V,. 

From this result, it follows that Hypothesis I is automatically satisfied, and thus, super- 
fluous whenever stated. 

Proof. First, the claims about ru(A) are obvious from the definition of along with 
Propositions ^, 12 and 19. 



From Proposition 19 together with equation ( |4.3| ), it follows that exp^2 A = ^2 (A) = 1 
for A £ T>c- Proposition |l^ implies that, for A £ int(/C), we obtain exp^2 A < 1) ^'^d thus 
by Proposition |ll|, r/(A) < 1. 



Also, Proposition 13 provides that, for each j £ A, the functions Xj — > ^2-\ ~^ 
r/(A) are continuous for A £ int(/C). Hence for d > 1, because Vc C IC, the function Xj r/(A) 
is left-continuous for X £ IC, and so is ^2 A- other words, r/(A) < 1 for A £ int(/C) and 
r/(-) is continuous and strictly increasing for A £ int(/C) with ?7(A) | 1 as A G /C approaches 
X £ Vc ioT d > 1. For d = 1, due to monotonicity, r]{X) = 1 for A G /C^, thus, t/(-) is constant 
and certainly right-continuous. The statement about the discontinuity is an instance of 
Proposition The statements pertaining to ^2;A carry over to Pj{X) and bj{X) for j £ A. 
Similar reasoning leads to the results for "ifp-x for each p > ru{X). □ 

Note that, due to the fact that ^2('^) = 1 for A G T^c, it follows that ru(A) > 2 for A G 
thus, in view of Corollary || and the left-continuity, that ru(A) = 2 for X £ T>c and d > 1. 

5.6 Weak Survival Region 

Corollary 6 For A G TZinlZ^, we have OiiX) < 1. 

Proof. For d = 1 indeed, by Corollary ^, 6*1 (A) = 02(A) = 1 for A G as desired. 

Next consider the case d > 1. Note that TZinTZI C JC because, for A G Vc, each bi{X) < 1. 
Hence by Corollary |5|, 0i{X) is left-continuous in each variable Xj for A G TZi D 1Z\. Now 
Lemma proves our claim for d> 1. □ 

Corollary 7 The weak survival region 7^2 has nonempty interior for d> 1. 

Proof. This is an immediate consequence of the following facts. Observe as in the previous 
proof that n ^ C /C. By Corollary |, 6*1 (A) < 1 for A G n By Corollary |, 
e2{X) = lioi X£Vc = V'c. In view of Proposition |l|, for A G n 

02(A) < 0i(A) < 1. 
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Consequently, 'Dcri{TZimZi) = 0. Furthermore, because 02(A) is continuous in each variable 
for A G int(/C), left-continuous for A G /C, and strictly increasing along directions of increase 
on /C, we conclude that the complement of TZi in fC has nonempty interior. 
Combining these with Proposition |l^ gives 

and Dc = ^3 H TZ^. This implies that /C \ TZi = IZ2, with IZ2 having nonempty interior. □ 
The discussion of the last proof has the following corollary. 

Corollary 8 For d > 1, we have 

/C^ = 7^3. (5.31) 



It remains to see what happens for d = \. Combining Proposition 19 with the left- 
continuity and strict monotonicity of /3(A) implies that (3 strictly increases in a continuous 
fashion all the way up to take the value 1 at I?c- Hence, the function /3(A) is continuous on 
(0, co). Since the contact process does not survive as long as /? < 1, it follows that IZ^ C IC^. 

We return to more than one dimension. For d > 1, the complement of /C coincides 
with the strong survival region TZ^ and T>c coincides with IZ2 H IZ^. As a consequence, 
7^2 n 7^3 C 7^.2, which leads to the following observation. 

Corollary 9 The symmetric anisotropic contact process on for d > 1 with A G 7^2 H 7^3 
survives weakly. 

This proves Theorem |2|. Replacing /C by TZi U 7^2 in all our previous results completes 
the proof of Theorem ^ and combining with Corollary ^ accomplishes Theorem 



Corollary 10 The statements of Proposition li hold for each A G 7^iU7^2 with each A^ > 



and the statements in Propositions and It hold for each A G intiJZi U7^2)- Furthermore 



for d > I, if continuity is replaced by left- continuity, then the statements in Propositions IS 



and It hold for each A G TZi U 7^2 with each > 0. 



Proof. It remains to see the last claim. But this is a consequence of the fact that i]{X) < 1 
for A 7^.3, and thus, r/(-) is continuous and strictly increasing for A G mt{TZi U 7^2) with 
r/(A) t 1 as A G 7^.2 approaches 7^.2 n 7^3 for d > 1. Hence, for each j G A, the function 
Xj is left-continuous for A G 7^2- But this implies left-continuity for A G 7^2 for all 

functions that are continuous for A G int(7^i U 7^2)- D 

For d = 1, the functions described in Corollary |l^ are constant on 7^3, thus certainly 
right-continuous for A G 7^3. 

Corollary 11 For d > 1, we have 

TuiX) = 1 for X £TZinTZ2, 

r„(A) = 2 /or AG ^2 n ^3- 
Moreover, for every X G TZi D TZ2, 

0i{X) = 1. (5.32) 
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Proof. Indeed, we have already noticed that ru{X) = 2 for A G Dc = 'R'2 H 7^3 for d > 1 
(remark following Corollary ||). As observed earlier, r„(A) > 1 for A TZi. By Lemma 
we have ru{X) < 1 for A G TZi n 7^2- All claims now follow from continuity and strict 
monotonicity considerations. □ 

Next we complete the proof of Theorem ^. 

Proposition 20 The weak survival region 7^2 enjoys the following properties: 
(a) The boundary TZi H TZ2 consists of all A such that the 6j(A) satisfy 

and the boundary 7^2 H TZ^ consists of all A such that the bi{X) satisfy 



1 



(b) Every line in the interior of the first quadrant in R*^ that passes through the origin 
has an intersection with IZ2 that is a line segment. 

(c) The region IZ2 has positive d- dimensional Lebesgue measure. 

(d) The region IZ2 is connected and is a symmetric region in the d parameters Aj. 



Proof. First we list some observations. Let d > 1. Define the functions R{x) = x/(l + x) 
and Qp{\) = J2ieA'^(^iW^)- '^'^^ function R{x) is continuously differentiable and has 
first derivative 1/(1 + x)^ > for nonnegative x, thus, is strictly increasing in x. Recall 
from (|]3|) that 6*1 (A) = 1 for A G ^1 n "^2 and from Proposition || that 612(A) = 1 for 
A G ^2 n'^s. 



(a) . Whence by (gj), for every A G 7^l n 7^2, 

Qi(A) = l, 

and for every A G 7^2 H 7^3, 

Q2(A) = 1. 

(b) and (c). In fact, (c) will follow once (b) has been proved. Fix A G TZi n 7^2- Thus 
by Proposition [l^, (32(A) < 1. For each a = {ak}keA with at least one ak > 0, due to the 
continuity properties of the functions bi{-), there exists a smallest real > such that 

Q2{{Xj + OjiJie^) = 1- 

The positivity of now guarantees claims (b) and (c). 

(d). The first portion of (d) is an immediate consequence of the statement in (b), 
whereas the second portion follows from the symmetric role that is being played by each 
function R(bi{X)^), i G ^, in defining Qi and Q2 together with the statement in (a). □ 
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The question whether, for d > 1, the critical contact process at the first phase transition 
behaves as in the lower or upper phase is more subtle than the analogous question at the 
second phase transition. At the second phase transition, the discontinuity settles the issue 
since the phases are well separated in a certain sense, whereas at the first phase transition 
the contact process behaves continuously. However, we exploit the fact that for X ^ TZi, the 
expected value of -E|3^m| must grow without bound as m — > oo, and thus, a Galton- Watson 
tree may be embedded in the set of vertices ever to be infected that is dominated by the 
contact process and whose attached Galton- Watson process is supercritical. 

The next result will prove Theorem |3|. 

Proposition 21 TZi n 7^2 C TZi, that is, for X £ TZi H 7^2, the contact process on for 
d > 1 almost surely becomes extinct. 



Proof. Let d > 1. Assume that Aa > A;, > Aj for all i £ A. Thus, A^ > and /?b(A) > 0. In 
view of (^), Lemma |5| and Corollary ^, we obtain J2x<=c%^ Wx ^ oo as m —>■ oo if and only 
if ExGS^ ^ oo as m ^ oo. 

It suffices to show that the Galton- Watson trees fr with mean offspring numbers J^xeOr '"a;' 
constructed in parallel with the Galton- Watson trees Tr, have corresponding subcritical 
Galton- Watson processes {Zn{r)}n>o for X £ TZi n7^2- Indeed, observe that the Galton- 
Watson process associated with dominates the set of infected vertices, in particular, on 
the event of survival, the limit set of fr contains the limit set A of the contact process. 
Thus, if the Galton- Watson process dies out, then the limit set A will be empty and the 
contact process cannot survive. 

The infection is most likely to follow trails along vertices with addresses composed of 
mostly letters a and b. For every integer m, for each real < s < 1, define Hg C Gm to be 
the set of all vertices at distance m from the root vertex which contain sm letters b and 
(1 — s)m letters a. We first determine the cardinality |-ff<i|. Indeed, note that \Hs\ = (^^) 
and expand by means of the Stirling formula for factorials, namely, ml = [m / e)^ ^ 2'Kme^ 
with (12m + 2/5m)"i < c < (12m)"^ to get 

\Hs\ = 




^27rs(l - s)m ^s*(l - s) 



Thus, by the mean value theorem along with the fact that each Ux{X) decays exponentially 
in |x| because A ^ 7^3, there is some < u^:{X, s) < 1 such that 



Ux{X) = -^^^^=^^^= 
xGHs V2vrs(l - s)m 



u^{X, s) 



(5.33) 



Set a(A,s) = n,(A, s)/(s^(l - s)!"^). 

Now let A G TZi n 7^2- If J2x€g,n ^xi^) — > oo as m — > oo (so that m is running through 
powers of 2), then it must be the case that there exists some s > such that J2xeHi, '^xW 
oo as m ^ oo. In that case, we must have a(A, s) > 1. By ( ^.21| ), the continuity and strict 
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monotonicity properties discussed earlier, there would exist some A' G int(7^i) in a direction 
of decrease for A such that 

1 < a(A', s) < a{X, s). 

Hence, it would follow that J2x£gm "^xi^') > Z^xeHs ^a;(A') ^ oo as m — > oo, equivalently 
by Proposition Yl,xi^c* '^xi^') — > oo as m — > oo. By Lemma ^, for sufficiently large 
m > 0, the Galton- Watson process Zn{m) attached to the Galton- Watson tree Tm would 
then be supercritical, and thus, the Galton- Watson tree Tm would be infinite with positive 
probability. Since this Galton- Watson tree is dominated by the set of vertices ever to be 
infected, this contradicts A' € TZi. Hence, we conclude that a(A, s) < 1 for all < s < 1. 
Finally, since a(A, s) < 1, an elementary calculation together with ( p. 33 ) yields that, for 



all sufficiently large m, the sum J2x£Hs '^x ^ I \/'^t^s(\ — s)m, which ^ as m ^ oo. 
Therefore, for each fixed 1 > s > and all sufficiently large m, each sub-Galton- Watson 
tree, restricted to the vertices in Hs, has attached a subcritical Galton- Watson process, thus, 
survives with probability zero. This implies that the full Galton- Watson tree f^, which is 
dominating the contact process, has an associated sitfecri^ica/ Galton- Watson process Z„(m). 
Consequently, the contact process cannot survive with positive probability. This completes 
the proof. □ 

Finally, recall that fin{V) = P{{i^n,'-^n+i, • • •) ^ ^} every ui = LO1LO2 . . . £ 0, and any 
Borel set y C and that the function 99 : S — > R was defined by ip{. . . X1X2 ■ ■ ■) = log 6x1 • 
Then the following result is an immediate consequence of Proposition ^ and Corollary |^. 

Theorem 10 For A € 7^2 cLnd every n > 1, the measure fin is absolutely continuous with 
respect to /z<^ and ^ l^tp as n 00. Furthermore, the pressure P{ip) = 0. 

6 Upper Bounds 

Now we are ready to prove the easier bounds for the Hausdorff dimensions. Recall from the 
beginning paragraph in Section 4 that 9i = 9 is the lead eigenvalue of Mi. 

Upper bound for the HausdorfT dimension 6{X) of the limit set. 
Lemma 10 With probability 1 on the event of survival, 6{X) < — log ^i/ log a. 



Proof. Recall that 3^m. denotes the set of all vertices x G Qm that are ever infected. Pick 
e > 0. Then by Lemma |5|, for all sufficiently large integers m, 

E\ym\ = J2 

x£g,n 

< 0i(A)(i+=)™. 

It follows from the Chebyshev-Markov inequality and the Borel-Cantelli lemma that almost 
surely, eventually 

\ym\<{ei{x)'+' + er. (6.1) 



57 



Observe that the sets provide a sequence of open covers of A, in particular, if 
denotes the set of all ends of T that pass through x, then 



For each x G y^, the diameter of £x (in the metric) is a"*. Hence, by (6.1) for all 
sufficiently large m, 

^ diameter„(f,)-'°s(^i(^)'+^+2-)/'°s- < 1 

(see e.g. Q). Because e > was arbitrary, this implies that with probability 1 on the event 
of survival, the Hausdorff dimension of A is 

o(A) < , 

log a 

as required. □ 

Upper bound for 5(A; /i). Let /i be an ergodic, cj-invariant probability measure on the 
space of semi-infinite reduced words from A. Recall that 5{\] fi) denotes the Hausdorff 
dimension of An (in the metric da), where is the subset of 17 consisting of all sequences 
u> that are "generic" for /x in the sense of definition ( 1.10| ). Recall that c^;^ : ^ R is the 
function defined by fx{xiX2 ■ ■ ■) = log6xi(A). Since ipx is continuous on relation ( 1.10 ) 
holds with f = ip. 

Lemma 11 For every e > 0, there exist sets Tm = Tmi^J-) C Qm of vertices at distance m 
from the root 1 such that 

lim -logir^l < h{fi)+e, (6.2) 

I m _ 

lim sup I — V log 6a;, (A) - / (pxd^i\<e, (6.3) 

and such that for every sequence X1X2 . . . G ^2, 

X1X2 . . . £ X1X2 ■ ■ ■ Xm £ rm(/i) infinitely often. (6.4) 

Proof. The Shannon-McMillan theorem and the ergodic theorem along with the ergodicity 
of /i guarantee that there are sets Tm C Qm such that 

lim /i( I J Vlix)) = 1, (6.5) 
lim max I — log ^(r2(x))/m — /i(/i) I = 0, (6-6) 



and such that ( |6.3| ) holds. A consequence of ( |6.6| ) is that for every e > 0, for sufficiently 
large m, and every x G Fm, 

exp{-m(l +e)/i(^)} < /i(f7(x)) < exp{-m(l - e)/i(^)}, (6.7) 

thus, 

|r,„|exp{-m(l+e)/i(^)} < ^( |J Q(x)) < |r„| exp{-m(l - e)/i(^)}, (6.8) 

X&Tm 
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which together with ( |6.5| ) provides (|6.2|). It remains to verify that ( |6.4| ) then holds. Pick a 
subsequence {mn}n>o of integers such that 

oo 

U n{x))]<^. 

n=l xeFmn 

In view of the Borel-Cantelli lemma, for ^u-almost every sequence x = xiX2 ■ ■ ■ G f^, for 
sufficiently large n, the initial segment X1X2 ■ ■ ■ Xm„ is contained in r^^. This completes the 
proof. □ 

Our next result proves half of (|1.12| ). 

Proposition 22 // + / ipxdfi < 0, then with probability one, A D = 0. // 
^(a*) + / 'fxdfJ' > 0, then with probability one on the event of survival, 

5(A;^)<-M±il^. (6.9) 
log a 



Proof. Let = TmifJ-) be as stated in Lemma 11, If a reduced semi-infinite word 
cj = X1X2 ... is an element of A n Q^, it must be the case that X1X2 ■ ■ ■ Xm G 3^m H Tm for 
infinitely many integers m. Define 

^m(^) = W = X1X2 . . . e Ann^ : X1X2 . . . a^m G 3^m n Tm}- 
Consequently for each m > 1, 

AnQf^C U Anifi). (6.10) 

n>m, 

Hence, the set Un>m^n(^) is a covering of A„(/i) by sets of diameter a™. Therefore, 
in order to find an upper bound for the Hausdorff dimension of A„(//), it suffices to find 
an upper bound for the cardinality of ym H Tm- By the same reasoning as in the proof of 



Lemma 10, for every e > and sufficiently large m, 



E\ymnTm\<{i+er E n^w = (i+^r E e^p{E'^A(^^x)}. 

XlX2. . .Xm&Tm j = '^ xeTm i = l 

From inequality (6^), it follows that |Fm| < exp{m(/i(;u) + e)}. Furthermore, by ( |6.3D , 



for every word x G F^, Z]j=i < m{J (pxdfi + e). With this in mind, we conclude 

that, for all sufficiently large m, the expected cardinality of ym H F^ is no greater than 
exp{m(/i(/i) + / ipxdfi + Se)}. The Borel-Cantelli lemma thus implies that, with probability 
one, eventually 

\ym n F^l < exp{m(/i(^) + J ipxdn + 4e)}. (6.11) 

In case h{fj.) + / ipxd^ + ^e < 0, eventually, ym<^Tm is empty. Thus, Am(/i) must be empty, 
and so by ( |6.10D , A n 0^ = 0. In the other case when h{fi) + / y^xd^ + 4e > 0, then as a 
consequence of inequality ( |6.11 ), for every n> 1 and all sufficiently large m, the set A„(^) 



is covered by exp{m(/i(/i) + / (pxdfi + 5e)} sets of diameter a"*. Since e > can be chosen 
arbitrarily small, 

SniKifJ-)) < -(/i(^) + J ^xdfi)/loga. 
Since A n 0^ C Un>m^n(M)i tii6 required inequality (|6.9|) follows. □ 
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7 Lower Bounds 



To verify formulae ( |1.7D and ( 1.12| ) for the Hausdorff dimensions of the random sets A 



and A n 0^, we need to estabhsh the lower bounds for the Hausdorff dimensions. For this 
purpose, we shall again consider Galton- Watson trees, embedded in the set of vertices of 
T that are ever infected, whose limit sets are subsets contained in A and A n 17^ that 
have Hausdorff dimensions which approach the required bounds. In turn, these Hausdorff 
dimensions are calculated by invoking a theorem of Hawkes 0, 21 ] and an extension given 



in ||16[. First, we show that the Hausdorff dimensions of both limit sets are almost surely 
constant. 

Lemma 12 The Hausdorff dimensions 5(A) and 6{X; ^) are almost surely constant on the 
event of survival. 



Proof. The reasoning is parallel for each of both Hausdorff dimensions. We shall give an 
argument for 5(A), largely borrowed from |15|. Let 5^, be the essential supremum of the 



random variable 5(A). Then for any 5 < 5,,, there is positive probability p that the limit set 
of a contact process initiated at the root 1 has Hausdorff dimension at least 5. Since the 
distribution of the Hausdorff dimension of a subset of the boundary dT does not change by 
an isometry of T, due to the geometry of the tree, it follows that for any vertex x G T, there 
is positive probability p that the limit set of a contact process initiated at x has Hausdorff 
dimension at least 5. Therefore, if J-t is the u-algebra generated by the percolation structure 
up to time t, then 

P{5{\)>5\Tt}>pI{\A,\>i}. 

But the martingale convergence theorem implies that this conditional probability converges 
to the indicator function of the event {5(A) > 5} almost surely as t ^ oo. Since, obviously, 
the indicator of the event {\At\ > 1} converges to that of the event that the contact process 
survives, it follows that 

/{5(A) >5} > P-f{survival} 

almost surely. Since the indicators take only the values and 1, it immediately follows that 
5(A) > 5 almost surely on the event of survival. Hence, 5(A) = 5* almost surely on the 
event of survival, as desired. □ 

7.1 Hawkes' Theorem and an Extension 

Recall the Galton- Watson trees r from Section 2.4 

Theorem 11 (Hawkes [0]) If the offspring distribution has mean /i > 1 and finite second 
moment, then almost surely on the event of nonextinction, the limit set Aqw of the Galton- 
Watson tree t has Hausdorff dimension (in metric da) 

dH[J^GW) = • 

log a 
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Hawkes examines only the case a = 1/2 but the result and its proof hold for all a G (0, 1) 
1 21, Now recall from Section 2.3 that the labelled Galton- Watson processes have label 



set B. For each label i £ B, let qi = J2FcB ieF Qi^) denote the probability that label i is 
included in a random set with distribution Q. Define a function tp : B^ ^ R by 

tp{xiX2...) =logqxi. (7-1) 

For any ergodic, shift-invariant probability measure // on the sequence space B^ , define Bj^ 
to be the set of //-generic sequences oj such that for every continuous function / : B^ —>■ R, 



1 " r 
lim - y /(a'co) = / fdfi. (7.2) 

1=1 



Theorem 12 (Lalley and Sellke [|16(] ) Letr he the labelled Galton-Watson tree attached 
to a supercritical labelled Galton- Watson process with label set B and offspring distribution 
Q, and let fj, be any ergodic, a-invariant probability measure on B^ . If h{fi) + J il'dfi < 
then with probability one, 

OtDB^ = 0. (7.3) 

If h{fi) + / ipdfi > 0, then almost surely on the event of nonextinction, the Hausdorff 
dimension of dr n B^ in the metric da is 

5^(9rn^,^) = -^^4±I^. (7.4) 
^ log a 

7.2 Lower bounds for 5{X) and S{X; jj) 

The following corollary together with Lemma ^ completes the proof of equation (1.7). 



Corollary 12 With probability 1 on the event of survival, (5(A) > — log (A)/ log a. 
Proof. The Hausdorff dimension 6h{^) of A is almost surely constant on the event of sur- 



vival by Lemma 12, thus, it is enough to show that for any real number 5* < — log0i(A)/a, 
there is positive probability that (5// (A) > 5* . Consider the embedded Galton-Watson trees 
Tfc. Since A € 7^2, we know that ^i(A) > 1. Indeed, if A is off the region TZ^ and in 7^2, 
then A € int('7^2), because TZi n 7^2 C TZi. Then the strict monotonicity of 9i{-) along any 
direction of increase and the fact that 6i{X) = 1 for A G T^i n TZ2 imply that Oi{X) > 1 
for A G 7^2- The second moment condition in the Hawkes' theorem is satisfied since the 
offspring random variable is a sum of indicator variables 

$Z ^{3 downward infection trail {root}— > x beginning at t=0} 

for every k. This sum of indicators obviously has finite variance for each k. By the Hawkes' 
theorem, on the set of nonextinction, the limit set A^ of has Hausdorff dimension (in the 
metric da) 

Sk = 

log a 
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Since the probability of nonextinction is positive, it follows that, with positive probability, 
the set A has a subset of dimension 5k- By ( ^.11 ), for sufficiently large k, 6k > S*, thus. 



with positive probability, the Hausdorff dimension of A is at least 5* , as desired. □ 
Proposition ^ together with the following result finishes the proof of equation ( |1.12| ) . 

Corollary 13 For any ergodic, shift-invariant probability measure ^ on $7, with probability 
one on the event of survival, 

<5(A;m) > -iKii) + I ipxdfi)/loga. (7.5) 



Proof. The proof is nearly parallel to the one of Corollary 12. Let ijj = (p in In the 



case that the righthand side of ( |7.5| ) is negative, the inequality apparently holds. Suppose 
otherwise. Consider the limit set A^, of the Galton- Watson tree r^. As in the previous proof, 
the underlying Galton- Watson process is supercritical and the second moment condition in 
the Extended Hawkes' theorem is satisfied. By the shift-invariance and ergodicity of /i and 
the definition of (p\ together with the Extended Hawkes' theorem, almost surely on the 
event of survival, the intersection of A^ with fi^ has Hausdorff dimension (in the metric da) 

h{fi) + /c-i ^iogn-=iCd/^(x) 



log a 

where the bj denote the entries of the matrix associated with the functions Wx (as 
opposed to Bi associated with the Ux, that we were using before). But by Proposition 
the matrices -Bf and Bi have the same lead eigenvalue, and thus, we conclude that bj = bj 
for every j £ A. Since C A, the last display is a lower bound for the Hausdorff dimension 

ofAnri^. □ 



8 Backscattering Inequalities 

Throughout this section, we only address the case d > 1. The term "backscattering in- 



equality" was coined in since the inequality 6{X) < 6h{^)/'^ may be proved by a 
backscattering argument in different context. The following argument is short and does not 
involve the backscattering idea. Another line of reasoning will be explained below that is 
based on the Gibbs Variational Principle. 

Proposition 23 For all X 7^3, 

OiiXf <92{X){2d-l) <2d-l, (8.1) 

with strict inequality Ol{X) < 2d — 1 holding except possibly for X G 7^2 H TZs. 
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Proof. Recall that 6i{X) and ^2 (A), respectively, are the lead eigenvalues of the matrices 
Mi(A) and M2(A), respectively. Thus, for i = 1,2, ^^(A) = lim^^oo(l*Mi(A)'"l)i/™, where 
1 denotes the 2d- vector with all entries 1. As an appeal to the Cauchy-Schwarz inequality 

(m 
ili2---im j=l 

(\ l/2m / \ l/2m 

E uhm El ■ 

The last line of the display tends to ^J62{\)^2d — 1 as m — > 00 because the number of 
reduced words iii2 ■ ■ - im of length m equals (2d — 1)™". Taking the limit as m — > 00 also 
implies that ^i(A)^ < 92{X){2d—l). It follows from Proposition |l^, Corollary and the strict 
monotonicity of 02(A) along any direction of increase that 02{X) < 1 with strict inequality 
except for A G 7^2 H TZ^. □ 

Corollary 14 For all A TZ^, with probability one, 

6iX) < hnm, (8.2) 
where strict inequality holds except possibly for A G 7^2 H TZ^. 



1 j jfi 



Proof. It is an easy exercise that the Hausdorff dimension of Q in the metric da equals 
— log(2d — l)/logQ!. Set 9 = 61. Hence, the result follows from equation (1.7) and Proposi- 
tion 12^. 



□ 



Gibbs Variational Principle. Consider any ergodic, shift-invariant probability mea- 
sure fi on 0, and the measure-preserving system {Q,a,fj,). Let h{fi) denote the measure- 
theoretic entropy, let 93 : — > R be any Holder continuous function, and let P{(p) denote 
the thermodynamic pressure of the potential function ip, as described in Section 3. The 
Gibbs Variational Principle |3|, ^ states that 

/i(/i) + / ^pdfi< P{ip), (8.3) 
Jn 

where the inequality is strict unless /i = /x<^ is the Gibbs state for if. If we let ip = p^px, 
where ip\{xiX2 ■ ■ ■) = log6i.^(A), the pressure functional is given by 

P(P(^a)= lim -logl fihiMY] =^oge{p-\). (8.4) 

\iii2...im j=l I 

Since the potential functions depend only on the first coordinate, the corresponding Gibbs 
states /Up^^ are the probability distributions of the stationary Markov chains with transition 
probabilities 

P(p;% = |^(i-5.(r^)), 
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where v is the lead right eigenvector of the Perron- Frobenius matrix Mp(A) and denotes 
the Kronecker delta function. Observe that these Gibbs states coincide with the ones 
discussed in Section 3. In particular, when p = 2 and A € 7^2 H TZs, this Gibbs state 
coincides with the measure announced in Theorem ^ (and discussed after that result). 
The following result now completes the proof of Theorem ^. 

Proposition 24 For every ergodic, a-invariant probability measure n on Q, 

6{\;ij) <^dH{%) (8.5) 

almost surely, where the inequality is strict unless A G 7^2 H TZs and p = p^. 



Proof. Let p = 2 and if = 2ipx. We combine (|8.3| ) and (|8.4D , divide both sides of the 



inequality by 2, and recall formula ( 1.12 ) for the Hausdorff dimension 6{X; p) of A n il^, 



h{p) + 2 f ifxdp < log ^2 (A) 
Jn 

^ h{p)+ ( ^xdp < -(/i(^) + log 02(A)) 
JQ. 2 



=> 6{X;p) < -^{h{p)+loge2{\))/loga. 

Unless p coincides with the Gibbs state P2if), = A**, strict inequality holds. Recall again 
that 02(A) < 1 for A in the interior of 7^2 and 02(A) = 1 for A € 7^2 H TZ^. This brings along 

<J(A;/x) < -^h{p)/loga, 

and strict inequality holds unless p = p^ and A G 7^2 H TZ^ . The statement now follows from 
the observation that fi^ has Hausdorff dimension, in the metric d^, equal —h{p)/loga, 
thanks to a theorem of Billingsley. □ 

Along with Corollary 14, the following proposition ends the proof of Theorem |5|. 



Proposition 25 Strict inequality holds in (l.i) for every A G 7^2 H 7^.3 except when the 



contact process is isotropic, that is, when all infection rates Aj = A. 

Proof. By Proposition it suffices to verify that 

0i(A)2 < 02(A)(2d-l) 

for every A G 7^2 H 7^3 except when the contact process is isotropic. Let A G 7^2 H 7^3. 
By applying formula ( |8.4D and the Gibbs Variational Principle to the potential functions 
1^ = ipx and if = 2ipx, respectively, with Gibbs state p = pi, we find by (^) 

2h{pi)+2 j ifxdpi = 21og0i(A) (8.6) 

h{pi)+2 j ^xdpi < log 02(A), (8.7) 
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with strict inequality in ( |8.7| ) except when the Gibbs states /xi and IJ.2 coincide. An instance 
of Theorem 1.28, allows this to happen when the difference of the potential functions (p\ 
and 2ipx is (co-) homologous to a constant function, equivalently, if and only if there exists 
constants e > and c > such that 

for every m > 1 and all finite reduced words iii2 ■ ■ - im £ Gm- This is possible if and only 



if all of the values bi, i A, are identical. However, by Proposition 18, this happens if and 
only if the contact process is isotropic or each bi = 1. But the latter would contradict that 
A G ^2 n TZs, thus, in 7^2. □ 

9 Isotropic Contact Process 

We first outline the existence proof for an immediate phase in the isotropic case prior 
to summarizing additional consequences in the isotropic setting. Indeed, we recall that 
A — > is a strictly increasing and continuous function as long as /3(A) < 1, left-continuous 
everywhere, that by Corollary /3(A) < l/{2d — 1) for A G TZ^ D TZi, and by Proposition 
19, /3(A) = l/\/2d — 1 for A S Vc- Thus, there must exist an intermediate interval between 
extinction and strong survival that has positive length. 

Corollary 15 For d > 1, consider the isotropic contact process, that is, with Aj = A, in the 
weak survival phase Ai < A < A2. Let /3(A) = /3i(A). Then there is some < Dx < 00 such 
that 

l/{2dDx,) < A < l/{2V2d^Dx,), 
l3{X) = ei{\)/{2d-l) and J(A) = - log[(2d - 1)/3(A)]/ log a, 

where 

0i(A) = l/{2\Dx) - {1/(4A2Z)2) _ {2d - 

In addition, 

^i(A) = (2d-l)/3(A), e2(A) = (2d-l)/3(A)2, ^2(A) = /3(A)ei(A). 

Especially at the phase transitions, 

OiiXi) = 1, /3(Ai) = l/(2(i-l), 

0i(A2) = \/2d^, /3(A2) = 1/V2d^, 
02(Ai) = l/V2d^, 

02{X2) = 1. 



Proof. Suppose that d > 1. (For d = 1, while equation ( 5.26D holds, relation ( 5.24| ) that 



will be relied on below is not valid since the set I?c is not in /C.) Note that continuity allows 



65 



to include the first critical value A = Ai in our calculations below. Also, observe that, in 
the isotropic case Aj = A, by Proposition 18, bi{X) = 6(A) = /?(A). 

First, it is an easy exercise to solve equation ( ^^ ) for (3, thus, P = 61/ {2d — 1). With 
this in mind, relation ( 5.24| ) 

l3 = \Dx[l + fP{2d-l)] (9.1) 

may be restated as 

ei/{2d-i) = \D^[i + el/{2d-i)i 

where there are some positive finite constants Ci and C2 such that for each Ai < A < A2, 
Ci < D\ = D\{x) < C2{n + 1) by Proposition 17. This quadratic equation in 9i has two 
solutions 

''"-' = 2Ak^\/5«! 

the relevant root being the one with the negative sign in front of the radical because /3 is a 
nondecreasing function in A, so is 9i. 

Next we find the endpoints of the interval (Ai,A2], the interval of the weak survival 
phase. Since ^i(Ai) = 1, by Corollary |l^, solving (|4.3|) yields /? = l/{2d — 1), and, solving 
(U) gives A = Ai = l/{2dDx) = l/{2dDx^). Moreover, since 6'2(A2) = 1 by Corollaries | 
and ^, solving equation ( |4.3| ) with p = 2 yields [3 = \j\/2d— 1, and, again solving ( p.l[ ) 
provides A = A2 = Xj {2^/2d - \Dx) = l/{2^/2d - IDx^). Additionally, from equation ( PD 
it follows that 62 = which equals f3'^{2d — 1). This finishes our proof. □ 

Lower bounds for each Xi^, k £ A^, for A G T^i n 7^2 or A G 7^2 H TZ^ may be found by 
determining the critical values for the anisotropic branching random walk (for a description 
of an algorithm, see Section 3.2), whose population dominates the one of the contact 
process. 



10 Growth Profile 

This section will be devoted to the proof of Theorem ^. Recall that rt and Rt denote the 
smallest and largest distances among the infected sites x £ At and Nn{ns) denotes the 
number of vertices x G Ans at distance n from the root that are infected at time ns. 



Proof of Theorem ^. Recall from ( p. 41 ) that g<I>2,s = n(^ip2 = logC for some s > 0. 



For each p > 'r„(A), the analogous identity is ric^^s^p^s = T^C^p- Therefore, it follows from the 
continuity and monotonicity properties of $i,s;. and Tp^., in A and the fact that Tp-^.-^ = for 
each A G TZi H 7^2 , by Theorem ||, along with ( |3.4C1| ) that there exists at least one solution 
s of = 0. 

Proof of ( |1.19 ). Suppose that exp{$i^s} is the lead eigenvalue of the time-dependent 



equivalent of the matrix Mi defined via the entries of the matrix associated with $2,s ■ Fix 
s > so that (^i s > 0. Then by the same arguments that led to the result in Lemma ||, we 
find, by relying on the time-dependent functionals, that 

lim ( V n,,„.)^/" = lim {ENn{ns)Y''' = exp¥i,, (10.1) 
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for every s > 0. Hence, for any e > and all sufficiently large n > 0, 

ENn{ns) < exp{n(¥i,s + e)}. 
Therefore, by the Borel-Cantelli lemma and the Markov inequality, almost surely, 

limsup — logA'^ri('^s) < log(exp{<I>i^s + e} + e), 

n— >oo n 

thus, since e > is arbitrary, 

limsup — logA^„(?is) < '^i.s- 

n—*oo n 

To prove the reverse direction, we grow some time-dependent labelled Galton- Watson 
trees embedded in the set of vertices that are infected at time ns. The construction is 



parallel to the one described in Proposition 14, from which we conclude that 



= '^i,s (10.2) 

for every s > 0. Along the same lines as before (Corollary ^), we also conclude that the 
mean offspring numbers for the Galton- Watson trees at time ns satisfy 

lim ui/™ = exp¥i,. (10.3) 

Since Nn{ns) dominates the corresponding Galton- Watson chain, it follows that, on the 
event of survival of the Galton- Watson process, 

lim inf — log A'^„(ns) > logexp <I>i ^ = <I>i ^ 

because by picking m suitably large, (a) the probability of the event of nonextinction of 
the Galton- Watson process is as close as desired to the one of the event of survival of the 
contact process, and (b) the expected number of offspring of the Galton- Watson process in 
Gn is as close as desired to the expected number of vertices in Qn that are infected at time 
ns. This finishes the proof of ( 1.19D . 



Proof of ( IiTTtD and ([TTT^ ). We will mimic the sketch of proof in [^] . Note that we 



have seen above that the equation = has at least one solution. Let si be the smallest 
solution and S2 the largest solution. Moreover, for any interval (a, 6), let Nt{a,b) denote 
the number of vertices x £ At with at < \x\ < bt. 

It suffices to show that, almost surely on the event of survival, for any 6 > 0, eventually 
Nt(t/si + 6t, oo) = and Nt{0,t/s2 — 5t) = 0. Once this is verified, it will then follow 
that, almost surely on the event of survival, limsup^^o^ Rt/t < si and liminf^^oo ^t/t > S2, 
which, in view of the part of the proof we have already seen, will imply (1.17) and ( |1.1^ ). 



Fix 5 > 0. A moment's thought shows that it is enough to consider integer times 
t. Observe that for each s < si, we have <I>i^s < 0. Pick e > small enough so that 
^i,s + e < with s = si/(l + 6si) < si. For sufficiently large fixed t, the probability that 
Nt{t/si + 5t, oo) > 0, by ( |10.lD , is no larger than 

E ^^'.t ^ E exp{n(¥i,, + e)} 

n>t/si+5t x£Q„ n>t/s\+5t 
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with s = Si/{l+5si) < Si. Since $i^s+e is negative, the sum is bounded above by cexp{—^t} 
for some positive constants 7 and c (depending only on 6). Since J2t^o^^P{~7i} < the 
Borel-CanteUi lemma provides that, almost surely, eventually Nt{t/si + St, 00) = 0. Since 
S > was arbitrary, this proves one half. But in fact, the proof that, almost surely, 
eventually Nt{0,t/s2 — 6t) = runs in parallel. This completes our proof. □ 
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